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Abstract 


The  focus  of  the  effort  being  carried  out  under  the  AFSOR  grant 
is  to  fully  characterize  the  behavior  of  rotating  stall  phenomenon  in 
turbomachinery  systems  and  to  assess  effectiveness  of  nonlinear 
control  strategies  using  both  CFD  simulations  and  experiments  on  the 
axial  compressor  facility  in  the  LICCHUS  lab.  Towards  accomplishing 
this  objective,  a  CFD  model  has  been  developed  and  it  has  been 
validated  using  experimental  data  available  in  the  literature.  The 
axial  compressor  facility  in  the  LICCHUS  lab  will  be  used  for 
correlation  with  CFD  simulations  and  for  stall  controller  studies. 


1  Introduction 


A  renewed  interest  in  modeling  and  simulation  of  aircraft  engine  components  has  emerged  during  the 
past  twenty  years.  One  motivation  for  such  interest  is  the  prediction  of  unsteady  stall  phenomena 
in  compressors  and  its  relation  to  fluid  dynamic  stability.  Modeling  of  complex  turbomachinery 
systems  has  grown  to  embrace  the  field  of  modern  control  theory,  in  search  for  mechanisms  of 
instability  suppression. 

The  pursue  of  simplified  models  stems  from  the  motivation  of  reducing  the  dimensionality  and 
associated  complexity  of  complete  aerodynamic  simulations.  While  these  efforts  offer  the  advantage 
of  computational  economy,  the  simplifications  made  result  in  the  residualization  of  physicaJ  system 
dynamics  which  might  be  of  interest,  and  whose  detailed  study  could  reveal  finer  details  of  insta¬ 
bility  inception  and  development.  However,  models  that  reduce  to  systems  of  ordinary  differential 
equations  possess  an#  additional  advantage.  The  impressive  advances  in  the  numerical  analysis, 
simulation,  and  control  of  complex,  ODE-based  nonlinear  systems  have  allowed  the  development 
of  control  strategies  for  turbomachinery  flow  instabilities  based  on  such  models.  In  this  manner, 
the  reduction  or  elimination  of  the  design  stall  margin  through  active  control  has  become  a  major 
objective  for  addressing  instability  phenomena  in  turbomachinery. 

In  principle,  the  practical  validation  and  implementation  of  these  control  strategies  will  enable 
the  design  and  manufacture  of  more  robust,  high-performance  machinery.  However,  as  pointed  out 
earlier,  most  existing  models  lack  quemtitative  description  capabilities  for  most  physical  systems. 
When  such  models  are  employed,  the  physical  system  under  consideration  must  be  tailored  to  the 
model  by  means  of  system  identification  tools  and  a  certain  degree  of  empiricism,  given  the  standing 
ideal-behavior  assumptions  present  in  these  models.  Such  detailed  numerical  analysis,  transient 
parameter  adjustment,  and  map  identification  are  required  when  the  control  law  to  be  implemented 
in  a  physical  system  is  derived  from  a  mathematical  model. 

While  no  simplified  model  can  adequately  capture  the  behavior  of  a  physical  system,  it  is 
believed  that  the  element  of  empiricism  can  be  removed  from  a  model  ibrmulation  by  studying 
detailed  physical  aspects  of  the  modeled  system.  For  example,  in  a  post-stall  model  validation 
program  carried  out  at  Georgia  Tech  it  was  found  that  the  addition  of  input  terms  to  mimic  the 
behavior  of  rotor  disturbances  results  in  some  improvements  in  the  dynamic  response  of  the  modeled 
system.  Simplified  models  available  in  the  present  day  literature  do  not  address  the  disturbances 
created  by  the  rotors  and  blade  row  interference  effects.  While  these  issues  can  be  studied  through 
experiment,  the  advances  in  instrumentation  and  measurement  in  turbomachinery  have  yet  to  insure 
an  increased  understanding  of  the  flow  pattern.  On  the  other  hand,  it  has  been  established  in  recent 
times  that  flow  field  information  obtained  from  the  numerical  solution  of  the  fundamental  equations 


of  fluid  motion  in  their  different  levels  of  simplification  has  become  a  profitable  means  of  obtaining 
and  analyzing  data  for  complex  turbomachinery  flows.  Given  that  the  behaviors  which  can  be 
analyzed  by  detailed  flow  simulations  are  quite  often  those  which  have  a  large  impact  on  the  response 
of  these  complex  flow  systems,  one  can  identify  an  incentive  for  considering  such  calculations,  as 
they  provide  information  not  completely  conveyed  by  simplified  models. 

In  this  spirit,  the  work  that  has  been  performed  under  the  grant  will  take  on  the  renewed 
objective  of  demonstrating  the  feasibility  of  simulating  stall  phenomena  in  turbomachinery  based 
on  the  numerical  solution  of  the  Navier-Stokes  equations.  The  simulation  of  stalled  flows  in  turbo¬ 
machinery  is  a  virgin  area  of  research;  it  is  expected  that  this  work  will  be  one  of  the  first  few  in  its 
class.  In  addition  to  these  efforts,  the  numerical  calculations  will  be  used  in  conjunction  with  the 
one  of  the  simplified-model-based  instability  control  concepts  presented  in  the  growing  literature  of 
stall  control  in  an  attempt  to  quantify  the  degree  of  success  of  such  strategies. 

The  major  contribution  of  this  work  to  the  existent  literature  is  expected  to  be  the  practical 
demonstration  of  the  possibility  of  increasing  current  knowledge  of  rotating  stall  instabilities  by  the 
application  of  well  established  tools  in  Computational  Fluid  Dynamics  (CFD).  It  is  also  expected 
that  exploring  and  analyzing  simple  simulation  results  will  generate  ideas  for  better  guidance  on 
post-stall  modeling  of  rotating  machinery. 

2  Progress  in  1995 

During  the  past  year,  validation  of  a  simplified  post-stall  model  against  experimental  data  from  an 
axial  flow  compressor  rig  at  Georgia  Tech  LICCHUS^  has  been  carried  out.  The  focus  has  been  the 
validation  of  nonlinear  behaviors  associated  with  bifurcations  in  the  system  dynamics,  as  captured  by 
the  model.  Steady-state  validation  has  been  carried  out  through  nonlinear  programming  techniques. 
These  optimal  results  work  quite  well  but  good  matching  over  some  flow/pressure  regimes  is  often 
hard  to  obtain;  however,  the  technique  can  be  used  directly  in  a  robust  control  design  methodology, 
in  order  to  account  for  uncertainty  in  the  model.  In  addition,  the  control  approach  presented 
in  the  grant-supported  publications  “A  Simplified  Approach  for  Control  of  Rotating  Stall  Pari  I: 
Theoretical  Development  and  Part  II:  Experimental  Results  (AIAA  Journal,  (11)6,  1195  —  1223, 
1995)  has  been  tested  in  the  model,  to  complement  the  experimental  and  analytical  results  presented 
in  those  papers.  The  results  obtained  for  the  validated  model  under  controlled  operation  (system 
gains  and  parameters)  will  be  used,  along  with  the  numerical  simulation,  to  establish  the  effect  of 
this  pzirticular  controller. 

^Acronym  for  Laboratory  for  Identification  and  Control  of  Complex,  Highly  Uncertain  Systems. 


In  order  to  perform  the  main  task  of  simulating  stalled  flows  in  a  compressor  rotor,  a  two- 
dimensional  multiblock  Navier-Stokes  based  flow  solver  called  MFOIL  and  developed  by  Prof.  L.N. 
Sankar  at  Georgia  Tech,  has  been  modified  for  the  application  of  solving  steady  and  unsteady  flows  in 
cascades.  During  the  past  year  and  after  approximately  nine  months  of  work,  confidence  in  the  code 
has  been  established  by  comparing  its  output  to  experimental  data,  both  in  steady  and  unsteady 
flow  scenarios.  The  numerical  approach  is  based  on  the  well-known  Alternating  Direction  Implicit 
(ADI)  scheme  employing  an  Approximate  Factorization  (AF)  on  the  linearization  of  the  equations. 

Typical  simulation  results  are  shown  in  Figures  1  and  2.  Figure  1  shows  a  comparison  of  the 
blade  loads  on  a  NACA  65-410  profile  cascade  at  zero  angle  of  attack  (ai  =  0),  30  degrees  flow 
off  the  axial  direction  (^i  =  30®),  cascade  solidity,  (t,  of  1.25,  free-stream  Mach  Moo  =  0.085,  and 
Reynolds  number  Re  =  245,000.  In  these  results,  loading  is  reported  in  terms  of  total-to-static 
pressure  difference, 


n  _  PO  Pturjact 

where  po  is  the  stagnation  pressure  of  the  upstream  flow.  For  incompressible  flow,  S  is  simply 
1  —  Cp,  where  Cp  is  the  standard  static-to-static  pressure  difference  coefficient.  The  experimental 
data  comes  from  the  NACA  Report  no.  1368  published  in  1958  by  J.  C.  Emery,  L.  J.  Herrig, 
J.  R.  Erwin  and  R.  Felix  under  the  title  ^Sysiemaiic  Two-Dimensional  Cascade  Tests  of  NACA  65- 
Series  Compressor  Blades  ai  Low  Speeds^ .  Other  cases  in  this  report  compared  well  with  the  output 
of  the  program.  Minor  discrepancies  can  be  attributed  to  the  application  of  boundary  conditions  for 
the  curvilinear  grid  employed  and  perhaps  an  unaccounted  small  bending  of  the  blades  under  load 
in  the  experiment. 

Figure  2  shows  a  comparison  of  the  amplitude,  |Cp,  |  and  phase  relative  to  input  of  the  lifting 
pressure  coefficient  Cp.  =  Cp.^^.^  “^C’pio^er  of  a  flat  plate  cascade  in  which  the  blades  undergo  simple- 
harmonic  pitching  about  their  leading  edges.  The  phase  angle  between  the  motion  of  neighboring 
blades,  otherwise  known  as  the  interblade  phase  angle  (IB PA),  is  zero  in  this  case,  as  well  as  the  mean 
angle  of  attack  and  blade  stagger  angle,  7  =  /?i  —  ai.  The  free  stream  Mach  number  is  0.5  and  the 
pitch  angle  amplitude,  ap  is  0.01  degrees.  The  reduced  frequency  of  oscillation  based  on  semichord 
is  ib  =  wcl{2Voo)  =  0.5,  where  w,  c  and  Voo  are  the  circular  oscillation  frequency,  blade  chord 
length,  and  free-stream  speed,  respectively.  The  results  are  compared  to  the  output  of  LINSUB,  a 
numerical  program  written  by  D.  S.  Whitehead  for  the  calculation  of  linearized  unsteady  subsonic 
flows  in  cascades  and  presented  in  volume  1  of  the  AGARD  Manual  on  Aeroelasticity  in  Axial 
Flow  Turbomachines  (AGARDograph  no.  298,  1987).  The  output  of  MFOIL  varies  considerably 
more  than  that  of  LINSUB  over  the  chord  length.  The  discrepancy  is  a  result  of  the  grid  resolution. 


coupled  with  a  time-varying  spatial  accuracy  at  the  surface  as  a  result  of  the  currently-implemented 
grid  deformation  technique  for  the  blade  motion.  Otherwise,  the  MFOIL  output  is  seen  in  good 
qualitative  agreement  with  that  of  LINSUB,  Although  this  issue  could  be  addressed  by  implementing 
a  grid-deformation  method  which  does  not  allow  the  grid  near  the  solid  surfaces  to  change  in  time, 
the  exploratory  research  to  be  carried  does  not  call  for  this  type  of  blade  surface  motions.  Such 
motions  have  only  been  coded  to  validate  the  cascade-modified  MFOIL  program  in  unsteady  flow 
cases. 


3  Plan  for  1996 

The  primary  objective  for  this  research  period  will  be  the  successful  simulation  of  rotating  stall  in  an 
isolated  rotor  based  on  the  numerical  solutions  obtained  through  MFOIL.  To  that  end,  it  is  expected 
that  the  first  five  months  of  the  year  will  allow  a  study  of  the  parameters  that  will  be  required  to 
obtain  a  propagating  stall  scenario  in  such  rotor. 

Once  the  the  occurrence  of  rotating  staJl  is  established,  the  issue  of  comparing  simulation 
outputs  with  simplified  model  outputs  will  be  undertaken.  According  to  simplified  models,  it  is 
necessary  to  account  for  the  fact  that  nonlinearity  has  a  strong  effect  on  the  behavior  of  the  math¬ 
ematical  formulations  which  model  these  systems.  In  particular,  the  presence  of  finite  domains  of 
attraction  of  certain  operating  conditions  will  be  investigated,  given  that  this  concept  has  allowed  the 
development  of  a  stall  controller  by  0.  O.  Badmus  at  Georgia  Tech  and  documented  in  his  Doctoral 
Thesis.  This  controller  will  be  implemented  in  the  numerical  code  to  assess  the  effectiveness  of  the 
control  strategy  and  to  point  out  the  feasibility  of  improving  modeling  through  CFD  simulations. 
This  part  of  the  work  is  expected  to  take  up  the  second  part  of  the  year  and  into  the  year  1997. 


Figure  1:  Comparison  of  Experimental  and  Numerical  Results  for  Cascade  Loading  on  NACA 
65-410  profile,  Moo  =  0.085,  Re  =  245,000,  Oi  =  0,  /?i  =  30°,  a  =  1.25. 


Figure  2:  Comparison  of  LINSUB  (“0”)  and  Numerical  Results  (“X”)  for  Lifting  Pressure  CoefR- 
dent  Cp„  on  Flat  Plate  Cascade  Oscillating  About  Leading  Edge,  k  =  0.5,  Moo  =  0.5,  IBP  A  =  0. 
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ABSTRACT 

For  the  control  of  combustion  instability  in  a  ramjet  engine  by  fuel  injection,  a  math¬ 
ematical  model  of  Fung,  Yang,  Sinha  and  Menon  reduces  to  a  forced  wave  equation  with 
delayed  boundary  control.  The  first  and  only  step  in  the  analysis  that  was  completed  dealt 
with  the  imforced  problem  which  exhibits  instabilities  even  for  small  delays.  The  main 
contribution  of  our  analysis  wels  to  relate  these  instabilities  to  the  instabilities  in  difference 
equations,  a  problem  that  is  better  understood  than  the  wave  equation. 
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Delayed  Feedback  Control  of  Hyperbolic  PDE  and  Combustion  Instabilities 

The  Physical  Problem.  Combustion  instability  in  a  ramjet  engine  involves  nonlinear 
interaction  among  acoustic  waves,  vortex  motion  and  imsteady  heat  release.  The  instability 
is  often  observed  as  a  large  amplitude  pressure  oscillation  in  the  low  frequency  range.  At 
some  criticzd  Hmit  of  the  sunplitude,  there  can  be  structural  damage  due  to  fatigue  or  there 
can  be  engine  unstart.  The  latter  occurs  when  the  shock  in  the  inlet  duct  is  expelled  to 
form  a  bow  shock  ahead  of  the  inlet  and  is  one  of  the  most  serious  problems  in  developing 
an  operational  ramjet  engine. 

The  results  of  experimental  ([1],  [2])  and  numerical  studies  ([3])  indicate  that  the 
pressure  fluctations  in  the  combustor  grow  to  large  amplitude  low  frequency  oscillations 
when  the  unsteady  heat  release  in  the  combustor  is  in  phase  with  the  pressure  fluctuation. 
These  studies  showed  also  that  a  large  scale  vortex/fleime  structure  propagates  through 
the  combustor  at  the  same  frequency.  This  suggests  some  nonlinear  coupling  between  the 
shear  flow,  the  pressure  oscillation  and  the  imsteady  heat  release. 

There  have  been  many  attempts  to  suppress  these  instabilities  using  both  passive 
controls  (geometric  modifications)  and  active  controls  (acoustic  feedback  control  and  sec- 
ondjiry  fuel  injection  control)  (see  [4]  for  a  survey).  Acoustic  controllers  appear  to  be 
imprEictical  due  to  the  hostile  environment  in  the  combustor  ^d  the  large  power  require¬ 
ments.  Active  control  techniques  using  secondary  fuel  injection  have  been  shown  to  be 
effective  in  certain  situations  and  also  have  the  possibility  of  increasing  the  net  thrust  of 
the  engine  (see  [5]),  and  therefore  seem  to  be  a  promising  approach  to  the  control  of  full 
scale  ramjet  engines. 

Practical  application  of  secondary  fuel  injection  control  in  a  full  scale  combustor  has 
not  been  demonstrated.  Part  of  the  difficulty  arises  from  the  fact  that,  in  a  high  flow 
rate  combustor,  the  combustion  is  irregular  and  multiple  frequencies  are  amplified.  In 
this  case,  the  dominant  instability  mode  can  change  with  the  operating  conditions.  Even 
when  control  is  applied  to  suppress  the  dominant  mode,  additional  new  frequencies  can  be 
excited  and  become  unstable  (see  [6], [7], [8]). 

A  Mathematical  Model.  Our  study  is  concerned  with  an  active  control  based  on  fuel 
injection  based  on  a  model  developed  in  [9],  [10],  [11].  This  approach  is  based  on  a  wave 
equation  which  describes  the  dynamic  behavior  of  nonlinear  oscillations  with  distributed 
feedback  actions.  The  resulting  boundary  value  problem  is  a  forced  linear  wave  with  the 
forces  including  all  influences  of  acoustic  motion,  mean  flow  and  combustion  response  under 
conditions  without  external  forcing  together  with  the  effects  of  the  control  inputs,  one  of 
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which  is  a  source  term. 

A  closed  loop  control  system  is  obtsiined  by  relating  the  source  term  to  the  mass  flow 
rate  of  the  injected  fuel  by  means  of  a  time  lag  theory  developed  in  [12].  This  term  is 
proportional  to  the  rate  of  change  of  the  injection  rate  of  the  control  fuel  at  a  delayed  time 
and  a  spatial  distribution  fimction  characterizing  the  fraction  of  the  fuel  element  burned 
at  a  fixed  position  with  a  time  delay  with  respect  to  the  moment  of  injection. 

The  system  is  implemented  by  first  determining  the  state  of  the  acoustic  field  by 
monitoring  the  instantaneous  pressure  signature  at  a  fixed  location,  for  example,  at  the  wall 
near  the  downstream  diffuser  ([13]).  This  can  be  interpreted  as  a  point  sensor,  located  at  a 
fixed  position  with  a  certain  amplication  factor  c.  Following  the  idea  in  [12],  a  Proportional- 
plus-integral  (PI)  controller  is  then  introduced  to  modulate  the  mass  injection  rate  of  the 
secondary  fuel.  Since  the  time  derivative  of  the  fuel  injection  rate  exerts  direct  influence  on 
the  acoustic  field,  the  PI  control  law  is  equivalent  to  a  Proportional-plus-Derivative  (PD) 
control  law.  Therefore,  the  acoustic  pressure  input  produced  by  the  combustion  of  the 
injected  fuel  is  given  as  a  linear  function  of  the  error  and  its  rate  of  change  at  a  delayed 
time,  with  proportionality  constants  Kp  and  A'd.  The  time  delay  between  the  sensor 
output  and  fuel  injecton  accounts  for  the  actual  time  required  for  data  acquisition,  signal 
processing  and  dynamic  response  of  the  fuel  injection  mechanism. 

The  parameters  Kp  and  Kd  are  control  parameters.  In  addition  to  these  parameters, 
there  is  the  freedom  in  the  selection  of  the  time  delays  and  the  spatial  distribution  of  the 
external  forcing  associated  with  the  burning  of  the  injected  fuel. 

In  [11],  the  authors  approximated  the  above  problem  by  using  a  finite  number  of  modes 
and  gave  a  procedure  to  relate  the  control  parameters  Kp,Kd  and  the  time  delays  r][  to 
the  elimination  of  instabilities.  Numerical  results  also  were  obtained  which  showed  that 
the  reduction  of  the  amplitude  of  the  lowest  mode  was  approximately  the  same  as  obtained 
by  the  experimental  results  in  [6].  On  the  other  hand,  it  was  also  noted  numerically  that 
another  low  frequency  was  excited.  They  were  unable  to  explain  the  source  of  this  mode 
due  to  the  fact  that  the  simulation  was  not  carried  out  long  enough  to  obtain  suflBcient 
data  to  spectrally  resolve  this  low  frequency. 

Even  though  the  above  model  applies  only  to  acoustic  effects  (ignoring  all  convective 
components),  there  are  still  many  unanswered  questions  which  appear  to  be  fundamental. 
This  is  particularly  true  with  respect  to  the  meinner  in  which  stabilization  of  some  modes 
seems  to  lead  to  instability  of  other  modes. 

The  effects  of  delays  in  boundary  control.  The  above  model  in  combustion  reduces 
to  a  discussion  of  the  dynamics  of  a  wave  equation  when  it  is  subjected  to  a  delayed  control 
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at  some  points  on  the  boundary.  It  has  been  known  for  some  time  that  the  introduction 
of  the  delay  can  lead  to  unexpected  behavior  in  the  dynamics.  For  example,  for  the 
linear  unforced  wave  equation  for  which  instantaneous  boimdary  control  yields  exponential 
stability,  it  is  known  that  arbitrarily  small  time  delays  can  destroy  stability  (see  [14]  [15] 

116J). 

The  fact  that  phenomenon  of  this  type  could  occur  had  been  observed  long  before  these 
particular  control  problems  were  considered  (see  [17],  [18],  [19],  [20],  [21],  [22]).  Although 
these  earlier  investigations  were  concerned  mainly  with  neutral  delay  differential  equations, 
a  simple  change  of  coordinates  shows  that  the  above  hyperbolic  PDE  is  equivalent  to  a 
neutral  delay  differential  equation.  The  instabilities  are  a  consequence  of  the  dynamics 
associated  with  a  difference  equation  arising  from  the  neutral  delay  differential  equation. 

The  main  accomplishment  so  far  in  this  research  has  been  to  make  this  identification 
clear.  The  problem  is  still  under  investigation. 
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Abstract 

Part  one  of  these  lectures  is  devoted  to  two  fixed  point  theorems  motivated  by 
dynamics  in  delay  equations:  an  aympototic  fixed  point  theorem  which  be 
applied  to  determine  w-periodic  solutions  of  an  w-peiiodically  forced  equation 
and  an  ejective  fixed  point  theorem  which  can  be  applied  to  the  determination 
of  nontrivial  periodic  solutions  of  autonomous  equations.  Part  two  is  devoted  to 
large  delays,  Hopf  bifurcations  and  fixed  points  of  maps.  Part  three  shows  that 
small  delays  can  destroy  stability  properties  in  delay  differential  equations  as  well 
as  boundary  control  of  partial  differential  equatons. 


Key  Words:  Delay  equations,  dynamical  systems,  asymptotic  fixed  point  theo¬ 
rems,  ejective  fixed  point  theorems,  periodic  solutions,  stability,  singular  pertiur- 
bations,  boundarj'  control. 

Fixed  Point  Theorems  h^otivated  by  Dynamics 

1.  Asymptotic  Fixed  Point  Theorems.  Let  us  begin  by  considering  the  ordinary 
differential  equation 

(1-1)  x  =  /(t,x) 

where  /  :  R  x  R"  ^  R”  is  a  smooth  function,  /(t,x)  =  f{t  +  u,x)  for  all  t,i.  It  is  of 
interest  for  such  an  equation  to  know  if  there  exists  a  periodic  solution  of  the  same  period 
as  the  vector  field;  namely,  an  cj-periodic  solution.  If  the  equation  arises  as  a  model  of 
a  physical  system,  then  it  is  expected  that  all  of  the  solutions  will  be  defined  for  <  >  0 
and  also  that  they  remain  bounded.  Is  this  enough  to  imply  that  there  is  an  w-periodic 
solution?  In  general,  this  is  not  the  case  and  determining  conditions  imder  which  this  is 
true  has  led  to  interesting  asymptotic  fixed  point  theorems.  In  this  section,  we  survey  some 
of  the  known  results  on  this  problem.  We  point  out  also  the  types  of  conditions  that  will 
imply  similar  results  for  problems  in  infinite  dimensions  and  that  will  have  applications  to 
functional  differential  equations  and  partial  differential  equations. 

*  Partially  supported  by  NSF,  DARPA,  AFOSR 
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The  determination  of  u-periodic  solutions  of  (1.1)  is  equivalent  to  finding  the  fixed 
points  of  the  Poincare  map  T,  where  Tzo  =  x(u>,xo)  and  x(t,xo)  is  the  solution  of  (1.1) 
satisfying  i(0, xq)  =  Xq.  The  following  result  is  due  to  Massera  (1950). 

Theorem  1.1. 

(i)  If  n  =  1  and  theie  is  a  solution  of  (1.1)  bounded  on  [0,  oo),  then  there  is  an 
uf-periodic  solution. 

(ii)  If  n  =  2,  all  solutions  of  (1.1)  are  defined  for  t  >  0  and  there  is  a  solution  of  (1.1) 
bounded  on  [0,oo),  then  there  is  an  uj-periodic  solution. 

(iii) If  n  =  2,  then  there  is  a  vector  field  f  such  that  there  is  a  solution  of  (1.1)  bounded 
on  [0,oo)  and  there  does  not  exist  an  uj-periodic  solution. 

(iv)  If  n  =  3,  then  there  is  a  vector  field  f  such  that  all  solutions  of  (1-1)  are  bounded 
on  [0,oo)  and  there  does  not  exist  an  u-periodic  solution. 

The  proof  of  (i)  is  very  easy  because  the  Poincare  map  is  monotone.  In  fact,  if  x(t,xo) 
is  bounded  on  [0, oo),  then  the  sequence  T^xq,  k  =  1,2,...,  is  bounded  and  monotone. 
Thus,  there  is  a  limit  x*  and  it  must  be  a  fixed  point  of  T.  The  proof  of  (ii)  is  more 
difficult  and  uses  a  fixed  point  theorem  of  Brouwer  which  relies  heavily  upon  the  fact  that 
the  dimension  of  the  space  is  two.  The  construction  of  counterexeimples  to  the  result  for 
n  =  2,  n  =  3  as  stated  in  (iii),  (iv)  are  difficult  (for  another  example  for  n  =  3  due  to  S.-N. 
Chow,  see  Yoshizawa  (1975)). 

In  the  case  where  the  vector  field  in  (1.1)  is  affine,  it  is  possible  to  show  that  the  exis¬ 
tence  of  a  bounded  solution  implies  the  existence  of  an  w-periodic- solution  (see  Yoshizawa 
(1975),  for  example).  This  result  also  is  true  for  affine  mappings  T  on  infinite  dimensional 
spaces.  To  state  the  result,  we  need  some  additional  notation  which  also  will  play  an  im¬ 
portant  role  in  our  subsequent  discussion  of  nonlinear  problems.  For  any  bounded  subset 
5  of  a  Banach  space  A ,  the  Kuratowski  measure  of  noncompactness  Oi(B)  is  defined  as 

o:(B)  =  inf{d  :  B  has  a  finite  cover  U  of  diaimeter  <  d}. 

A  map  r  ;  A'  -+  A*  is  said  to  be  condensing  if  a(T{B))  <  a{B)  for  any  boimded  set  B 
in  X  for  which  a{B)  >  0.  If,  in  addition,  there  exists  a  constant  li:,0  <  ib  <  1,  such  that 
a(T(B))  <  ko:(B)  for  any  boimded  set  B  in  X,  then  T  is  said  to  be  an  a- contraction. 
The  following  result  is  essentially  due  to  Chow  and  Hale  (1974).  We  give  the  proof  since 
it  is  very  simple  and  uses  the  fixed  point  theorem  of  Massatt  (1980)  (which  generalizes 
a  result  of  Darbo  (1955)  on  or-contracting  maps)  and  states  that  a^  condensing  map  of  a 
closed  bounded  convex  subset  of  a  Banach  space  into  itself  must  have  a  fixed  point. 

Xlieorem  1.2.  If  X  is  a  Banach  space,  z  ^  X  is  fixed,  L  is  a  linear  Q~condensing  map 
on  X ,  Tx  =  Lx  +  z  for  X  e  X  and  there  is  aii  xq  e  X  such  that  {T^xo,k  =  0, 1, . . .}  is 
bounded,  then  there  is  a  fixed  point  of  T  in  X. 
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Proof.  Let  D  be  the  convex  hull  of  the  set  {z"  =  T"zo,n  =  1,2,...}.  If  y  €  P,  then 
y  =  E,€  j0ix\  where  J  is  a  finite  subset  of  the  positive  integers,  >  0,  =  1.  Since 

Ty  =  Ly  +  2  =  I(Eiej^,z*)  +  {'Liej0i)z 

=  ^i^j0i{Lx'  +  2)  =  T.i^j0iTx^  €  P, 

we  have  T(D)  C  D.  Since  T  is  o-condensing,  there  is  a  fixed  point  of  T  in  D. 

Let  us  return  to  the  nonlinear  system  (1.1).  Levinson  (1944),  in  his  study  of  the 
periodically  forced  \'an  der  Pol  equation,  initiated  the  study  of  the  modem  theory  of 
dissipative  systems.  If  T  is  a  continuous  map  on  a  Banach  space  X,  we  say  that  T  is  point 
dissipative  if  there  is  a  bounded  set  B  in  X  if,  for  any  x  £  X,  there  is  an  no  =  no(z,  B)  such 
that  T”z  €  B  for  n  >  no.  If  X  is  finite  dimensional  and  T  is  point  dissipative,  Levinson 
(1944)  showed  that  the  map  T  has  a  maximal  compact  invariant  set  A  (that  is,  TA  =  A). 
He  then  showed  that  the  Poincare  map  for  the  periodically  forced  van  der  Pol  equation  was 
point  dissipative  and  used  the  above  property  of  existence  of  a  maximal  compact  invariant 
set  to  show  that  there  is  an  integer  k  such  that  there  is  a  periodic  solution  of  period  kui-, 
that  is,  a  subhamionic  of  order  k.  From  Theorem  l.l(ii),  since  n  =  2,  we  know  that  there 
must  be  a  periodic  solution  of  period  u.  Other  results  allow  us  to  prove  the  following 
rather  remarkable  result. 

Theorem  1.3.  If  the  Poincare  map  for  (1-1)  is  point  dissipative,  then  there  is  a  fixed 
point  ofT. 

A  proof  of  this  result  can  be  supplied  (see,  for  example,  Pliss  (1966),  Yoshizawa  (1975)) 
using  the  following  asymptotic  fixed  point  theorem  of  Browder  (1959): 

Theorem  1.4.  Let  S  and  Si  be  open  convex  subsets  of  the  Banach  space  X  and  let  So 
be  a  closed  convex  subset  of  X  such  that  Sq  C  Si  C  S.  If  T  :  S  X  is  continuous,  T{S) 
is  contained  in  a  compact  set  of  X  and  if,  for  a  positive  integer  m,  T"*  is  defined  on  Si 
and  Uo<j<mT^'^o  C  Si  and  T'^Si  C  So,  then  T  has  a  fixed  point  in  So- 

It  was  also  observed  by  Gerstein  and  Kiasnoselskii  (1968),  Horn  (1970),  Billotti  and 
LaSalle  (1971)  that  Theorem  1.4  could  be  used  to  extend  Theorem  1.3  to  the  case  where 
T  is  completely  continuous.  Jones  (1965)  and  Yoshizawa  (1966)  used  a  similar  method  to 
obtain  u>-periodic  solutions  of  retarded  functional  differential  equations  (RFDE)  when  the 
period  to  is  at  least  as  large  as  the  delay.  This  condition  was  imposed  to  have  the  Poincare 
map  completely  continuous. 

Theorem  1.3,  when  extended  to  completely  continuous  maps  on  a  Banach  space,  is 
very  useful  in  the  applications  to  parabolic  partial  differential  equations  (parabolic  PDE) 
and  to  RFDE  when  the  period  u;  >  r,  where  r  is  the  delay.  What  happens  when  a;  <  r? 
What  about  other  applications  that  deal  with  neutral  functional  differential  equations 
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(NFDE)  and  damped  hyperbolic  PDE?  Extensions  have  been  made  which  are  sufficiently 
general  to  apply  to  these  situations.  To  describe  one  of  the  main  results,  we  need  some 
additional  concepts.  We  say  that  the  continuous  map  T  :  X  X  is  compact  dissipative  if 
there  is  a  bounded  set  B  in  X  such  that,  for  einy  comi^act  set  K  C  X,  there  is  an  integer 
no  =  no(A,B)  such  that  T”/\  C  B  for  n  >  no .  The  following  result  was  independently 
proved  by  Hale  and  Lopes  (1973),  Nussbaum  (1972). 

Theorem  1.5.  Jf  X  is  a  Banach  space  and  T  :  X  X  is  condensing  and  compact 
dissipative,  then  T  has  a  fixed  point. 

A  set  A  is  smd  to  attract  points  {attract  compact  sets)  (attract  bounded  sets)  of  a  set 
B  if,  for  any  c  >  0,  and  any  point  x  e  B  (compact  set  H  in  B)  (bounded  set  U  in  B), 
there  is  an  integer  no  such  that  T^x  (T"B)  (T’^U)  belongs  to  the  e-neighborhood  of  K  for 
n  >  Uq.  It  is  known  that  the  conclusion  of  Theorem  1.5  is  still  valid  if  we  replace  ‘compeict 
dissipative’  by  ‘there  is  a  compact  set  which  attracts  compact  sets  of  X'  (see,  for  example, 
Hale  and  Lopes  (1973),  Nussbaum  (1972),  or  Hale  and  Verduyn-Lunel  (1993)). 

We  give  a  brief  indication  of  the  proof  of  Theorem  1.5  following  the  ideas  in  Hale  and 
Lopes  (1973)  since  it  can  be  deduced  from  the  following  result  of  Horn  (1970),  which  is  an 
extension  of  Browder’s  Theorem  1.4. 

Theorem  1.6.  Let  So  C  Si  C  S^  be  convex  subsets  of  a  Banach  space  X  with  So,S2 
compact  and  Si  open  in  S2.  Let  T  :  S2  —*  X  be  a  continuous  mapping  such  that,  for  some 
integer  m  >  0,  PSi  C  S2  forO  <j  <m-l  and  T^Si  C  So  for  m  <  j  <  2m  -  1.  Then  T 
has  a  Bxed  point. 

To  use  this  result,  Hale  and  Lopes  (1973)  first  prove  the  following  intriguing 

Lemma  1.1.  Suppose  that  K  C  B  CZ  S  are  convex  subsets  of  a  Banach  space  X  with 
K  compact,  S  closed  and  bounded,  and  B  open  in  S.  If  T  :  S  X  is  continuous, 
{T^B,j  ^  0}  C  5  and  K  attracts  points  of  B,  then  there  is  a  closed,  bounded,  convex 
subset  A  of  S  such  that 


A  =  cb[Uj>ir-'(B  n  A)]  A  n  A"  7^  0.' 

Using  Lemma  1.1  and  Theorem  1,6,  it  is  possible  to  prove 

Lemma  1.2.  Suppose  that  K  C  B  C  S  are  convex  subsets  of  a  Banach  space  X  with 
K  compact,  S  closed  and  bounded,  and  B  open  in  S.  If  T  :  S  X  is  continuous, 
{T^B,j  >  0}  C  5,  AT  attracts  compact  sets  of  B,  and  the  set  A  of  Lemma  1.1  is  compact, 
then  there  is  a  Bxed  point  of  T. 

It  is  now  easy  to  see  that  T  being  a-condensing  implies  that  the  set  A  of  Lemma 
1.1  is  compact.  To  verify  that  the  conditions  of  Lemma  1.2  are  satisfied  if  T  is  compact 
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dissipative  requires  some  effort,  and  the  verification  uses  arguments  that  are  more  closely 
related  to  stability  theory  (see  Hale  and  Verduyn-Lunel  (1993)). 

As  a  simple  application  of  Theorem  1.5,  let  us  consider  the  RFDE 

(1-2)  x{t)  =  f{t,x,) 

where  X({6)  =  x(t  +  6),  0  €  [— r,  0],  r  >  0,  and  the  initial  data  is  chosen  to  be  in  the 
space  C  =  C([-r,0]  ;  R").  We  assume  that  the  function  f  :  B.  x  C  —*  IR"  is  smooth 
and  that  f(t  +  lo, ip)  =  f{t,(p)  for  all  t^ip.  The  Poincare  map  T  :  C  C  is  defined  as 
[Tip]{6)  =  x{t  +  0,V’)‘>  ^  G  where  x{t^ip)  is  the  solution  of  (1.2)  with  initial  data  ip 

on  [-r,0].  Let  us  assume  that,  for  each  fixed  t  >  0,  the  set  {x(s,v3),0  <  s  <  t,v?  G  B}  is 
bounded  if  B  is  bounded.  Since  a  solution  of  (1.2)  is  continuously  differentiable  for  t  >  0, 
if  u;  >  r,  then  the  Ascoli-Arzela  theorem  implies  that  T  is  completely  continuous.  If  a;  <  r, 
this  is  no  longer  true.  On  the  other  hand,  if  we  write  the  solution  z(-, y?)  of  (1-2)  as 

x<(-,V>)  =  S{t)p  +  U{t)ip, 

where  t/j  =  S{t)p  satisfies 

y(t)  =  0  for  t  >  0, 
y{t)  =  <p{t),  for  t  e  [-r,0], 

then  it  is  easv  to  see  that  U{t)  is  a  completely  continuous  operator.  Furthermore,  for  any 
7  >  0,  there  is  a  ^  >  0  such  that 

li5(0IU(C;C)  <  t  >  0. 

Therefore,  if  we  change  the  norm  in  C,  the  Poincare  map  T  is  an  o-contraction.  Also,  it  is 
not  too  difficult  to  show  that,  if  T  is  point  dissipative,  then  it  is  compact  dissipative.  We 
can  now  use  Theorem  1.5  to  prove  that  there  is  a  fixed  point  of  T  without  any  restriction 
on  the  period  uj. 

For  other  applications  of  Theorem  1.5  to  NFDE,  hyperbolic  PDE,  the  beam  equation, 
etc.  see,  for  example.  Hale  (1988). 

We  end  this  section  with  some  remarks  about  autonomous  systems;  that  is,  evolution¬ 
ary  equations  on  a  Banach  space  A'  which  generate  a  C° -semigroup  T{t\t  >  0.  In  this 
case,  for  any  r  >  0,  the  map  =  T{t)  is  a  Poincare  map.  If  we  let  Fix  (2V)  be  the  set  of 
fixed  points  of  and  we  know  that  Urg(o,i]  Fix(rr)  belong  to  a  compact  set,  then  it  is 
not  difficult  to  show  that  there  is  an  equilibrium  point  for  the  semigroup;  that  is,  a  point 
xo  €  X  such  that  T(t)xo  =  xo  for  all  t. 

To  state  a  precise  result,  we  use  the  concept  of  a  global  attractor.  A  set  in  AT  is  the 
global  attractor  for  T{t)  if  A  is  compact,  invariant  and  attracts  bounded  sets  of  X.  The 
positive  orbit  7'^(B)  of  a  bounded  set  B  in  X  is  defined  as  the  set  {r(t)x  :  t  >  0,x  €  B}. 
The  following  result  is  contained  in  Hale  (1988). 
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Theorem  1.7.  Suppose  that  T(t)  =  S(t)  +  I7(t)  :  X  -*  X  is  a  C® -semigroup,  C/(t) 
is  completely  continuous  and  there  is  a  continuous  function  k  :  IR'*'  — »  11"*^  such  that 
k(t,r)  — ♦  0  as  <  — »  oo  and  ||5(<)x||  <  k(t^r)  if  ||x|l  <  r.  If  T{t)  is  point  dissipative  and 
is  bounded  if  B  is  bounded,  then  the  global  attractor  A  exists  and  there  is  an 
equilibrium  point  ofT(t). 

We  only  remark  about  the  equilibrium  point.  The  conditions  imply  that  we  can  change 
the  norm  in  X  so  that  Tr  is  an  o-contraction  eind  then  Fix(Tr)  C  A  and  we  can  use  the 
remark  preceding  the  theorem. 

For  autonomous  ordinary  differential  equations,  the  existence  of  a  zero  of  vector  field 
can  be  proved  imder  weaker  conditions  as  in  (i)  and  (ii)  of  Theorem  1.1. 

-4n  interesting  example  due  to  Jones  and  Yorke  (1969)  for  autonomous  ordinary  dif¬ 
ferential  equations  in  shows  that  there  is  a  vector  field  for  which  all  solutions  are 
bounded  and  yet  there  is  no  zero  of  the  vector  field.  The  construction  makes  use  of  the 
fact  that  there  are  minimal  sets  in  whicli  do  not  separate  R^  into  two  open  sets  with 
one  bounded  and  the  other  unbounded  (in  particular,  a  torus  can  be  a  minimal  set). 

2.  Ejective  Fixed  Point  Theorems.  Let  us  first  define  what  we  mean  by  an  ejective 
point  of  a  map. 

Definition  2.1.  Suppose  that  is  a  Banach  space,  [7  is  a  subset  of  X  and  x  is  a  given 
point  in  U.  Given  a  map  T  :  U  \  {x}  — >  X,  the  point  x  e  U  is  saud  to  be  an  ejective  point 
of  T  if  there  is  an  open  neighborhood  G  C  X  of  x  such  that,  for  every  y  e  G  HU,  y  ^  x, 
there  is  an  integer  m  =  m(y)  such  that  T'^y  ^  G  DU. 

To  motivate  a  general  ejective  fixed  point  theorem,  we  consider  first  some  simple 
examples.  Consider  the  ODE 

(2.1)  i  =  /(i). 

where  x  G  R^  and  /  is  a  smooth  function  with  /(O)  =  0.  Suppose  that  it  is  known  that 
each  solution  osdllates  about  the  origin;  that  is,  if  x  =  (xi,X2)  and  x(i)  is  a  solution  of 
(2.1)  with  initial  data  xq  on  the  positive  xi-axis  K  =  {(ii.O)  :  Xi  €  [0,oo)},  then  there  is 
a  first  time  t  for  which  x(r)  €  K.  The  set  K  is  a  cone  and,  if  we  define  Txo  =  x(t),  then 
T  :  K  -*  K.  The  map  T  has  a  fixed  point  on  the  boundary  dK  of  K  given  by  x  =  0,  which 
corresponds  to  the  equiUbrium  point  0  of  (2.1).  A  nontrivial  fixed  point  of  T  corresponds 
to  a  periodic  orbit  of  (2.1).  Without  some  additional  conditions  on  the  map  T,  zero  may 
be  the  only  fixed  point.  If  we  suppose  that  the  origin  is  unstable,  then  there  exists  an  r  >  0 
such  that  |rx|  >  |x|  for  0  <  [x]  <  r;  that  is,  0  is  an  ejective  fixed  point  of  T.  Therefore,  if, 
for  0  <  |xj  <  r,  we  know  that  the  set  {T^x,  fc  =  0, 1,2, — }  is  bounded,  then  we  will  have 
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a  nontrivial  fixed  point  of  T.  This  will  be  true,  for  example,  if  we  suppose  that  \Tx\  <  |i| 
if  |xl  >  for  R  sufficiently  large.  This  is  a  simple  ejective  fixed  point  theorem. 

These  ideas  are  easily  extended  in  principle  to  (2.1)  where  x  €  IR^.  In  fact,  sup¬ 
pose  that  the  equilibrium  point  0  is  hyperbolic  with  one  real  negative  eigenvalue  and  two 
complex  eigenvalues  with  positive  real  parts.  Then  the  linearization  of  (2.1)  about  zero 
has  a  one  dimensional  stable  manifold  l'F*(0)  and  a  two  dimensional  imstable  manifold 
Ty“(0).  Let  us  choose  the  coordinates  x  =  (xi,X2,X3)  so  that  Vr'(O)  =  {(0,0,X3),X3  6  IR) 
and  H^“(0)  =  {(ii,X2,0),(xi,i2)  €  IR^}.  From  the  fact  that  the  origin  is  a  saddle 
point,  there  is  a  cone  K  with  nonempty  interior  (and  therefore  of  three  dimensions) 
which  contains  H^“(0)  \  {0}  and  is  positively  invariant  under  the  flow  for  the  lineariza¬ 
tion  near  0.  Also,  becaues  the  eigenvalues  with  positive  real  parts  are  complex,  the  set 
K  =  {x  €  A'  :  I  =  (xi,0,0),xi  >  0}  has  the  property  that,  for  any  xq  €  A",  there  is  a  first 
return  point  ToXo  to  A  under  the  linear  flow.  The  set  K  is  a  cone  with  0  being  an  ejective 
fixed  point  of  Tq. 

If  we  now  siq)pose  that  the  cone  K  also  is  positively  invariant  tmder  the  flow  defined 
by  (2.1)  and,  for  any  xq  €  A,  there  is  a  first  return  time  Txo  to  K  under  the  flow 
defined  by  (2.1),  then  we  have  T  ;  A  — >  A  and  0  is  an  ejective  fixed  point  of  T.  If  we 
suppose  also  that  there  is  a  A  >  0  such  that  |rx|  <  jx]  for  |x|  >  R,  then  we  can  define 
A/j  =  {x  €  A  :  |x|  <  A}  and  have  T  :  Kji  — ^  Kfi  with  0  being  an  ejective  fixed  point. 
Using  a  small  amoimt  of  theory  of  Liapunov  functions,  it  is  possible  to  deduce  from  the 
instability  of  zero  that  there  is  a  r  >  0  such  that  the  set  A^r  =  {a:  G  Ar  :  |x|  >  r)  has 
the  property  that  T  :  Kur  — »  A^r-  Brouwer’s  fixed  point  theorem  imlies  that  there  is  a 
fixed  point  of  T  in  A^r  and,  thus,  a  nontrivial  periodic  solution  of  (2.1). 

The  method  described  in  the  previous  paragraphs  actually  has  been  applied  to  a 
specific  problem  in  circuit  theory  concerned  with  the  quenching  of  undesirable  oscillations 
(see  Oldenburger  and  Boyer  (1961),  Hale  (1963)). 

It  is  natural  to  attempt  to  carry  the  above  procedure  over  to  infinite  dimensional 
systems  by  replacing  the  Brouwer  theorem  by  the  Schauder  theorem  for  compact  maps. 
Unfortunately,  it  is  not  so  easy  to  exclude  a  small  neighborhood  of  the  origin  using  Li- 
apimov  theory.  Convexity  (or  being  equivalent  to  a  convex  set)  is  not  easy  to  verify. 
Motivated  by  an  example  to  be  discussed  below,  Browder  (1965)  gave  an  ejective  fixed 
point  theorem  for  compact  maps,  which  was  later  extended  by  Nussbaum  (1974)  to  the 
following  result. 

Theorem  2.1.  If  A  is  a  closed,  bounded,  convex,  infinite  dimensional  set  in  X,  T  : 

\  {^o}  A  is  an  a-contraction,  and  xq  is  an  ejective  point  ofT,  then  there  is  a  fixed 
point  ofT  in  A  \  {xo}.  If  K  is  finite  dimensional  and  xq  is  an  extreme  point  of  A,  then 
the  same  conslusion  holds. 

By  combinipg  in  a  very  interesting  way  the  concept  of  eigenvalue  for  cone  maps  as 
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in  the  work  of  Krasnoselskii  (1964)  and  Grcifton  (1969),  Nussbaum  (1974)  proved  the 
following  more  global  result.  For  any  M  >  0,  let  Sm  =  {x  6  X  :  |x|  =  M},  Bm  =  {x  6 
X:|x!<M}. 

Theorem  2.2.  If  K  is  a  closed  convex  set  in  X,  T  :  K  \  {0}  — ♦  K  is  an  a -contraction, 
0  €  K  is  an  elective  point  of  K,  and  there  is  an  M  >  0  such  that  Tx  —  Ax,  x  €  /v  D  S^, 
implies  A  <  1,  then  T  has  a  hxed  point  in  K  D  \  {0}  if  either  K  is  infinite  dimensioned 
or  X  is  finite  dimensional  and  0  is  an  extreme  point  of  K 

In  specific  application  of  these  results  to  fxmctional  differential  equations,  it  is  possible 
to  take  advantage  of  some  of  the  geometric  properties  of  the  flow  in  order  to  obtain  the 
ejectivity  of  the  special  point  in  the  statement  of  the  above  theorems.  We  state  this 
precisely  in  the  proper  context. 

Let  us  now  turn  to  a  specific  delay  differential  equation,  referred  to  as  Wright’s  equa¬ 
tion  (Wright  (1955),  (1961)),  which  was  studied  in  detail  by  Jones  (1962)  and  was  the 
motivation  for  the  development  of  the  above  general  ejective  fixed  point  theorems. 
Suppose  that  q  is  a  positive  constzuit  and  consider  the  equation 

(2.2)  i(t)  =  -ax{t  -  1)[1  x(t)]. 

If  C  =  C([— 1,0],IR),  then,  for  any  6  C,  there  is  a  unique  solution  x(t,<,5,o)  of  (2.2) 
satisfying  x(^, (,?, a)  =  <p{6).,6  6  [—1,0].  If  we  assume  that  v?(0)  >  —1,  then  it  is  easy  to 
see  that  x(t,(,?,a)  exists  for  all  t  >  0  and  x(t,  97,0)  >  —1  for  all  t  >  0.  In  the  following, 
we  will  always  assume  that  <,9(0)  >  -1  and  still  use  the  notation  <p  e  C.  If  we  define 
[^o(0v^](^)  =  ^  €  [—1,0],  then  Ta(t),i  >  0,  is  a  C°-semigroup  on  C. 

Our  objective  is  to  determing  periodic  orbits  of  (2.2).  Let  us  try  to  mimic  the  proce¬ 
dure  outlined  above  for  (2.1)  for  n  =  3.  To  do  this,  we  first  must  understand  the  behavior 
of  the  solutions  of  the  linear  vziriational  equation  about  the  origin: 

(2.3)  x{t)  =  —oix(i  —  1). 

The  eigenvalues  of  this  equation  are  those  values  of  A  for  which  there  is  a  solution  e'^';  that 
is,  those  A  which  satisfy  the  characteristic  equation 

(2.4)  A  +  Qe-'^  =  0. 

This  equation  has  infinitely  many  solutions  with  Re  A  00  as  |A|  -♦  00.  Other  properties 
of  these  eigenvalues  is  given  in  the  following  result. 

Lemma  2.1. 

(i)  If  0  <  Q  <  7r/2,  every  solution  of  (2.4)  has  negative  real  part. 
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(ii)  If  o  >  c~*,  there  is  a  solution  A(o)  =  7(0)  +  *<7(Q')  of  (2.4)  which  is  continuous 
together  with  its  first  derivative  in  a,  7'(a)  >  0,  o\oc)  >  0,  0  <  a(a)  <  tt,  a(7r/2)  =  7r/2, 
7(7r/2)  =  0,  and  7(0)  >  0  for  a  >  ■jr/2. 

Lemma  2.1  implies  that  there  is  a  Hopf  bifurcation  at  o  =  7r/2.  It  is  possible  to  show 
that  it  is  supercritical  (that  is,  in  a  neighborhood  of  x  =  0,  a  =  7r/2,  there  is  a  periodic 
orbit  only  if  o  >  7r/2,  it  is  unique,  hyperbolic  and  stable)  (see  Chow  and  Mallet-Paret 
(1977)).  Since  aiv/'l)  =  7r/2,  the  bifurcating  periodic  solution  has  period  approximately  4 
and  is  slowly  oscillating-,  that  is,  the  distance  between  zeros  of  the  solution  is  greater  than 
the  delay  (>  1). 

As  a  increases,  there  are  more  eigenvalues  which  cross  the  imaginary  axis  and  lead 
to  Hopf  bifurcations.  However,  these  periodic  solutions  are  not  slowly  oscillating.  This 
suggests  that,  if  we  attempt  to  carry  out  a  program  similEir  to  the  mentioned  one  for  ODE, 
we  must  restrict  attention  to  the  set  SO  of  slowly  oscillation  solutions  of  (2.2). 

For  all  values  of  a,  the  dominant  imstable  eigenvalues  are  the  ones  in  part  (ii)  of 
Lemma  2.1.  The  corresponding  eigenfunctions  span  a  two  dimensional  subspace  W  of  C 
and  there  is  a  cone  K  in  C  with  nonempty'  interior,  W  \  {0}  belongs  to  the  interior  of  K 
and  K  is  positively  invariant  under  the  flow  defined  by  the  linear  equation  (2.3).  Since  the 
dominant  eigenvalues  are  complex,  solutions  of  (2.3)  rotate  around  K.  If  we  knew  that  the 
same  properties  hold  for  the  nonlinear  equation  (2.2),  then  we  would  be  close  to  situation 
in  ODE.  All  of  these  remarks  can  be  made  more  precise  and  we  do  this  following  the  ideas 
in  Jones  (1962).  However,  we  first  state  the  main  result  for  (2.2). 

Theorem  2.3.  If  a  >  n/'2,  then  (2.2)  has  a  nontrivial  periodic  solution  which  is  slowly 
oscillating. 

We  only  give  the  main  steps  in  the  proof  and  refer  to  Hale  and  Verduyn-Lunel  (1993) 
for  the  details.  We  first  need  a  result  on  the  oscillatory  properties  of  the  solutions  of  (2.2). 

Lemma  2.2. 

(i)  If  v?(0)  >  —1  and  the  zeros  of  x(*,  ip,  a)  are  bounded,  then  x(t,  (p,a)  0  as  t  00. 

(ii)  lf<p(0)  >  —1,  then  x{t,p,a)  is  bounded.  Furthermore,  if  the  zeros  of  x{-,(p,a)  are 
unbounded,  then  any  maximum  of  x(t,<p,a),  t  >  0  is  less  that  e“  —  1. 

(Hi)  If  v’(O)  >  —  1  and  o;  >  1,  then  the  zeros  of  x(',(p,a)  are  unbounded. 

(iv)  If(p(6)  >  —1,  — 1  <  ^  <  0,  then  the  zeros  of  x(t,<p,a)  are  simple  and  the  distance 
from  a  zero  of  x(t,(p,a)  to  the  next  maximum  or  minimum  is  >  1. 

Let  K  be  the  class  of  all  monotone  increasing  functions  <p  £  C  such  that  (p(9)  >  0, 
— 1  <  ^  <  0,  ^(—1)  =  0.  Also,  suppose  that  0  €  if .  The  set  if  is  a  cone.  If  o  >  1,  9?  €  K, 
V?  ^  0,  let 

z(v?,a)  =  min{t :  x(t,y>,ot)  =  0,x(<,</3,a)  >  0}. 
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This  minimum  exists  from  Lemma  2.2,  Parts  (iii)  and  (iv).  Also  >  2.  Furthermore, 

Lemma  2.2(iv)  implies  that  x(<,(,?,a)  is  positive  and  increasing  on  (2(9,  o),  2((,5,a)  +  1]. 
As  a  consequence,  if  t{^p,q)  =  q)  +  l,  then  we  can  define  the  mapping  A  •.  K  —*  K  by 

A(o)0  =  0 

A{a)ip  =  ra(r(v?))((^),  v’  7^  0, 
where  Ta{i)  is  the  semigroup  defined  by  (2.2). 

Lemma  2.3.  The  map  t  :  K  \  {0}  x  (l,oo)  — +  (0, 00)  defined  by  t(9,  o)  =  2(<p^a)  +  1  is 
completely  continuous. 

From  Parts  (ii)  and  (iv)  of  Lemma  2.2,  it  follows  that  \A{a)<p\  <  e“  - 1  for  each  €  K 
and  A{a)  takes  any  bounded  set  B  in  K  \  {0}  into  the  set  {y?  G  C  :  |v?|  <  e"  -  1}.  Using 
Lemma  4.3  and  the  fact  that  A(0)  =  0,  it  is  possible  to  show  that  -4(a)  is  continuous  at  0 
and  that  A(a)  is  completely  continuous. 

If  we  let  /iTcr  ={<,?€  A  :  |v^|  <  e“  —  1},  then  A  :  Ka  is  completely  continuous 

with  the  fixed  point  0  an  extreme  point  of  /vq.  Also,  the  set  I\a  is  a  closed  bounded  convex 
subset  of  C.  If  the  point  0  is  ejective,  then  we  can  use  Theorem  2.1  to  complete  the  proof 
of  Theorem  2.3.  The  proof  of  this  is  a  little  technical,  but  follows  from  the  fact  that  the 
dominant  part  of  the  local  unstable  manifold  of  the  origin  is  tangent  at  the  origin  to  the 
two  dimensional  manifold  W  of  the  dominant  part  of  the  unstable  manifold  of  the  linear 
equation  (2.3).  We  do  not  give  the  details. 

If,  in  (2.2),  we  let  1  +  x(t)  =  then  we  obtain  the  equation 

(2.5)  i:(t)  =  af{y{t  -  1)), 

where  /(y)  =1-6’'  which  satisfies  the  property  that  /(y)  has  negative  feedback;  that  is, 
y/(y)  <  0  for  y  /  0  and  /'(O)  <  0.  For  equations  of  the  form  (2.5)  with  some  additional 
conditions  on.  the  function  /,  the  above  method  has  been  applied  to  obtain  the  existence  of 
nontrivial  periodic  solutions  (see  Hale  and  Verduyn-Lunel  (1993),  Supplementary  Remarks 
to  Chapter  11). 

Equations  (2.5)  and  more  general  ones  of  the  type 

(2.6)  x{t)  +  <Tx{t)  =  f{x{t  -  1)) 

occur  very  often  in  the  applications.  The  following  interesting  result  is  due  to  Hadeler  and 
Tomiuk  (1977). 
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Theorem  2.4.  Suppose  that  a  is  a  positive  constant  and  f  has  negative  feedback.  If 
there  is  an  interval  I  such  that  f(I)  C  I  and  the  ongin  is  linearly  unstable,  then  there  is 
a  nontrivial  periodic  solution  of  (2.6)  which  is  slowly  oscillating. 

The  above  method  will  not  yield  directly  a  proof  of  this  theorem  since  the  slowly 
oscillating  periodic  solution  need  not  be  monotone  increasing  of  an  interval  of  length  one 
of  the  form  [2,2  +  1],  where  2  is  a  zero  of  the  solution.  The  important  new  idea  is  to 
replace  the  cone  K  used  for  the  proof  of  Theorem  2.3  by  another  consisting  of  fimctions 
which  when  weighted  by  a  specific  t3rpe  of  exponential  fimction  are  monotone  increasing. 

Let  us  concfider  the  following  generalization  of  (2.6)  which  has  been  used  as  a  model 
for  the  treinsmission  of  light  through  a  ring  cavity  (see  Vallee,  Dubois,  Cote  and  Delisle 
(1987),  Vallee  and  Mamott  (1989)) 

(2.7)  (S„^  ^  +  1)  •  •  •  (^1  ^  +  IMO  =  /(V(<  -  1))’ 

where  S  =  (Sj, .  ..,6^)  €  (0,00)*".  Hale  and  Ivanov  (1993)  have  used  the  ideas  of  Hadeler 
and  Tomiuk  (1977)  as  well  as  the  method  above  to  obtain  the  following  result. 

Theorem  2.5.  Suppose  that  I  is  a  bounded  interval  such  that  f{I)  C  I,  f  has  negative 
feedback.  Then  there  is  a  Sq  >  0  such  that,  for  each  £  G  (0, 6o)"*>  equation  (2.7)  has  a 
slowly  oscillating  periodic  solution. 

Theorem  2.5  is  Theorem  2.4  if  m  =  0.  For  m  =  2,  an  der  Heiden  (1979)  has  obtained 
the  conclusion  in  Theorem  2.5  with  any  smallness  restrictions  on  c  but,  of  course,  with  the 
additional  condition  that  the  origin  is  linearly  unstable. 
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Part  2.  Large  Delays  and  Oscillations 


1.  Introduction.  It  is  a  well  known  fact  that  increasing  the  delay  in  a  retarded  delay 
differential  equation  (RDDE)  often  has  a  tendency  to  destabilize  the  motion.  Also,  the 
destabilization  frequently  occurs  through  a  Hopf  bifurcation;  that  is,  as  the  delay  is  in¬ 
creased,  a  periodic  orbit  bifurcates  from  an  equilibrium  point.  We  have  encountered  this 
phenomenon  in  our  study  of  Wright’s  equation  (2.2)  of  Parti,  Section  2.  By  a  reseeding  in 
time,  the  parameter  a  in  (2.2)  plays  the  same  role  as  a  delaj'.  We  now  want  to  discuss  sit¬ 
uations  where  the  destabilization  occurs  through  a  Hopf  bifurcation  ajid  try  to  understand 
the  nature  of  the  profile  of  the  resulting  periodic  orbits  as  the  delay  approaches  infinity. 

We  begin  with  a  simple  example  of  a  linear  system  to  see  how  the  eigenvalues  behave 
as  a  function  of  the  delay. 

Example  1.1.  For  a.b  constants,  consider  the  RDDE 
(1-1)  i{t)  = —ax{t)  -  bx(t  —  1). 

As  in  Part  1,  the  initial  data  for  solutions  is  taken  in  the  space  C  =  C({-1, 0];  IR).  Also, 
if  is  the  solution  of  (1.1)  with  initial  data  <p  6  C  at  t  =  0,  then  we  define  the  C°- 

semigroup  T{t)  :  C  C  by  the  relation  T{t)(p[6)  =  x(t  -1-  6^ip)^6  £  [—1,0].  The  operator 
T{t)  is  compact  for  t  ^  1,  the  spectrum  of  T{t)  consists  onlj'  of  point  spectrum  and  it  is 
given  by  where  A  is  a  solution  of  the  characteristic  equation 

(1-2)  A -1- a -F  6e"^  =  0. 

In  the  (a,  6)-plane,  the  stability  diagram  is  show  in  Figure  1.1. 


Figure  1.1.  Stability  region  for  x{t)  =  -ax{t)  -  bx{t  -  1). 

In  the  region  between  the  line  b  =  -a  and  the  curve  asymptotic  to  the  line  a  =  b 
which  contains  the  wedge  jaj  <  |6|,a  >  0,  each  solution  of  (1.2)  has  negative  real  part  and 
the  zero  solution  of  (1.1)  is  exponentially  stable.  On  the  line  b  =  —a,  {a,b)  ^  (—1,1), 
there  is  a  simple  zero  of  (1.2).  On  the  cur\-e  asymptotic  to  the  line  a  =  6,  (a,  6)  ^  (-1, 1), 
there  are  two  purely  imaginary  roots  ±ta;(a,6)  with  tu(a,6)  -♦  tt  as  (a,  6)  -+  (oo,oo).  At 
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the  point  (a,  6)  =  (  —  1,1),  there  is  a  double  zero  root  with  nonsimple  elementary  divisors 
(see,  for  example,  the  Appendix  of  Hale  and  Verduyn-Lunel  (1993)). 

For  T  >  0,a,6  constants,  consider  the  equation 

(1.3)  x{t)  =  —a.x{t)  —  bx{t  —  r), 

which  by  a  change  of  time  scale  is  equivalent  to  the  equation  (1.1)  with  a  =  dr,  6  =  6r. 
If  we  suppose  that  (d,  6)  is  in  the  stability  region  for  (1.1),  then  we  can  destabilize  the 
origin  by  increa.sing  the  delay  if  and  only  if  6  >  |d|.  The  destabilization  occurs  by  having 
two  eigenvalues  cross  transversally  the  imaginary  eixis  at  a  critical  value  T(d,6)  >  0.  If 
we  subject  (1.1)  to  nonlinear  perturbations  which  are  second  order  near  the  origin,  then 
the  Hopf  bifurcation  theorem  implies  there  must  be  periodic  orbits  in  a  neighborhood  of 
(x,r)  =  (0,  r(d,i).  b 


Figure  1.2.  Stabilit}'  of  the  origin  of  (1-3)  for  (a,  5)  =  (1, 1  +  A)). 

To  be  more  specific,  let  us  assume  that  d  =  1,  i  =  1  +  A,  where  A  is  a  real  psuameter. 
In  this  case,  we  can  destabilize  the  origin  by  increasing  t  if  A  >  0  and  cannot  if  A  <  0. 
.Also,  if  A  >  0  and  if  we  define  ^(A)  =  r(l,  1  +  A),  then  ■r(A)  — >  oo  as  A  — »  0  and  the  purely 
imaginary  solutions  ±tu;(A,  r(A))  of  (1.2)  aie  such  that  u;(A,t(A))  — ^  tt  as  A  -♦  0.  Note 
that  the  corresponding  eigenfunction  is  peidodic  of  period  winch  approaches  2  as 

A  ^  0. 

If  we  keep  d  =  l,  fc=l  +  A  and  formally  put  r  =  oo  in  the  rescaled  version  of  (1.3), 
we  obtain  the  mapping  on  R;  i  •-+  —(1  +  A)i.  The  fixed  point  zero  of  this  map  is  stable 
for  A  <  0  and  unstable  for  A  >  0.  Furthermore,  if  nonlinear  terms  were  included,  then 
there  will  be  a  period  doubling  at  A  =  0. 

Let  us  now  consider  the  nonlinear  equation 

(1.4)  ei(t)  +  i(t)  = /x(x(t  -  1)), 

where 

(1-5)  /a(x)  =  — (l  +  A)i  +  (?(x),  ^(i)  =  O(i^)  as  I -b  0. 
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Recall  that  e  =  i/r,  where  r  is  the  delay. 

Let  (A,  c( A))  be  the  curve  in  the  (A,  €)-plane  along  which  there  are  two  purely  imaginary 
solutions  of  (1.2)  and  the  remaining  ones  have  negative  real  parts.  For  a  fixed  small  value 
of  A.  the  point  c(A)  is  a  point  of  Hopf  bifurcation  with  respect  to  e  for  (1.4).  In  the 
(A.€)-plane.  the  Hopf  bifurcation  curve  divides  the  upper  half  plane  intersected  with  a 
neighborhood  of  the  origin  into  two  regions  S  and  U  with  the  property  that  the  origin 
is  hyperbolic  stable  in  S  and  hyperbolic  unstable  in  U  (see  Figure  1.3)  provided  that  we 
reatiict  our  discussion  to  initial  data  from  the  subspace  SO  of  C  which  corresponds  to 
slovh  oscillating  solutions.  Recall  that  a  slowl}'  oscillating  solution  is  one  for  which  the 
distance  between  zeros  is  at  least  as  large  as  the  delay,  whicli  in  our  equation  is  1.  In  the 
subspace  SO.  the  dimension  of  the  unstable  manifold  of  the  origin  in  U  is  2. 

lestnct  out  attention  to  SO  because,  for  a  fixed  value  of  A,  there  are  a  countable 
numbei  of  values  of  c  at  w-hich  two  eigenvalues  cross  the  imaginary  axis  from  left  to  right 


Figure  1.3.  The  stable  and  unstable  regions  of  the  origin. 

Let  X  be  a  slowd}-  oscillating  periodic  orbit  of  (1.4)  in  a  small  neighborhood  of  z  =  0 
lor  f.  A  small.  €  >  0.  The  problem  that  w'e  w-ant  to  discuss  is  the  follow'ing: 

For  a  fixed  value  of  X,  what  is  the  limiting  profile  of  this  solution  as  e  — ♦  0.^ 

From  our  remarks  above,  we  know  that  the  period  is  approximately  2  since  the  eigen¬ 
values  ±iu>(A,  r(A))  on  the  Hopf  cmwe  have  the  property  that  a;(A,T(A))  -♦  tt  as  A  — »  0. 

For  e  =  0  in  (1.4),  w^e  obtain  the  map  x  fx{x).  It  is  natural  to  conjecture  that 
there  should  be  some  relationsliip  between  the  dynamics  of  (1.4)  near  the  origin  and  the 

map  fx  near  the  origin.  In  particular,  we  might  expect  the  periodic  orbits  to  be  related  to 
period  tw'o  points  of  the  map. 

The  point  A  =  0  is  a  point  of  bifurcation  to  period  two  points  of  the  map  fx.  Let  us 
suppose  that,  in  a  small  neighborhood  of  z  =  0,  A  =  0,  there  are  only  a  finite  number  of 
period  two  points  (a),  fif),  qA  >  0,  j  =  1,2, . . . ,  all  hyperbolic,  and  ordered  as  a}  < 
iff  <  k.  Define  (o:^,^o  )  =  (OtO)-  The  period  two  points  of  fx  will  alternate  their  stability 
properties;  that  is,  if  is  stable,  then  and  (a^+„i9/+i)  are  unstable. 

To  each  (q;A,^A )^  we  can  define  a  square  wave  2-periodic  function  Sj(t)  by  the  relation 
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^>(0  -  Qj  (resp.  fif)  if  0  <  <  <  1  (resp.  1  <  <  <  2)  (see  Figure  1.4a).  This  function 
IS  stable  (resp.  unstable)  for  the  map  fx  if  the  period  two  point  is  stable  (resp. 

unstable)  for  the  map  fx. 

When  the  period  two  point  (o^,^^),  j  >  is  unstable,  the  square  wave  2-periodic 
function  sj  has  an  infinite  dimensional  unstable  manifold.  Therefore,  we  do  not  expect  it 
to  be  related  to  any  periodic  orbit  that  begins  near  a  Hopf  curve.  Because  of  the  nature  of 
the  solutions  of  the  characteristic  equation  (1.2),  there  is  a  two  dimensional  center  manifold 
of(1.4)  for  each  (A,c)  near  (Ao,c(Ao)).  If(A,e)  is  small,  we  also  expect  that  there  should  be 
some  two  dhnensional  invaiiant  manifold  that  contains  all  slowdy  oscillating  periodic  orbits 
near  the  ongm  with  each  periodic  orbit  on  this  manifold  encircling  the  origin.  If  each  such 
periodic  orbit  is  hyperbolic,  then  they  alternate  between  being  stable  and  unstable;  that 
IS,  the  index  of  the  corresponding  Poincaie  map  is  either  zero  or  1.  It  is  to  be  expected 
under  reasonable  conditions  that  the  limiting  periodic  orbits  as  «  0  should  be  related  to 

the  period  two  points  of  the  map  fx. 


Figure  1.4a.  Square  wave  sj  Figure  1.4b.  Pulse  wave  pj 

.A.S  a  consequence  of  the  above  intuitive  remarks,  when  the  period  two  point 
j  >  1,  IS  unstable,  we  define  a  pulse  2-periodic  function  p^(t)  bv  the  relation  p^(t)  = 
(resp.  D})  if  ,  =  0  (resp.  (  =  1).  p>(,)  =  if  0  <  (  <  1  (resp' 

1  <  t  <  2)  (see  Figure  1.4b).  The  pulse  2-periodic  function  is  unstable  for  the  map  fx  and 
the  dimension  of  the  unstable  manifold  is  1  (if  we  exclude  translations  in  t  by  a  constant). 
Following  our  intuition  that,  for  c  small,  we  expect  that  the  orbits  of  (1.4)  near  the 

origin  should  behave  in  some  way  as  the  orbits  of  the  map  />,'we  make  the  following 
conjecture,  ^ 

Conjecture.  For  e,  A  small  and  in  a  sufficiently  small  neighborhood  of  the  origin,  for  each 
J  >  1,  there  exists  a  periodic  solution  if  of  (1.4)  of  period  approximately  2  such  that 

(1)  -  sJ  0  as  e  -V  0  uniformly  on  compact  subsets  of  IR  \  Z  if  (o^,;3^)  is  stable 

(2)  xf  -  pA  0  as  c  -»  0  uniformly  on  compact  subsets  of  El\  Z  if{Q^,^\)  is  unstable, 

wfiereZ=  {n:n  =  0,±l,±2,...}.  ^ 

If  the  intuitive  remarks  made  before  the  statement  of  the  conjecture  could  be  made 
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precise,  one  could  probably  prove  the  conjecture.  Unfortunately,  we  do  not  know  at  this 
time  how  to  do  this. 

.A.nother  possible  approach  to  a  proof  could  be  the  following.  At  each  point  (Aq,  €(Ao)) 
on  the  Hopf  bifurcation  curve,  it  is  possible  to  use  the  method  of  Liapunov-Schmidt  to  ob¬ 
tain  the  bifurcation  function  G(A,  c,  c)  for  (A,  c)  in  a  neighborhood  of  (Ao,  €(Ao))  and  for  the 
approximate  amplitude  c  of  the  periodic  solutions  of  period  approximately  27r/u;(Ao,c(Ao)), 
where  ia,’(Ao,  e(Ao))  is  the  pureh'  imaginary  eigenvalue  on  the  Hopf  bifurcation  curve.  We 
know  that  a»(Ao,e(Ao))  — »  s-  as  Aq  — »  0.  The  zeros  of  G(A, e, c)  correspond  to  periodic 
solutions  of  (1.4)  with  period  close  to  2.  Therefore,  it  is  to  be  expected  that  these  zeros 
should  be  related  to  the  periodic  points  of  period  2  of  the  map  fx;  that  is,  to  the  fixed 
points  of  the  map  f^.  If  this  could  be  proved,  we  would  be  very  close  to  a  proof  of  the 
conjecture. 

For  the  case  in  which  there  is  a  generic  period  doubling  for  the  map  fx,  the  above 
conjecture  is  known  to  be  true.  Let  us  state  the  result  more  precisely.  Suppose  that 

(1-®)  /a(2')  =  -(1  +  A).r  -1-  ax'  -|-  bx^,  0  =  +  b  ^  0. 

We  lemark  that  it  is  not  necessai'v  to  assume  that  fx  is  a  polj'nomial  in  x  or  that  the 
dependence  upon  A  is  as  simple  as  in  (1.6).  In  fact,  it  is  enough  to  assume  that  fx  has  the 
property  that 

fx{0)  =  0  for  all  A, 

5x/o(0)  =  -l,  dlxfo{0)  =  -l, 

dlfoiO)  =  a,  ^/o(0)  =  b. 

We  assume  (1.6)  only  for  some  simplicity  in  notation. 

Under  the  assumption  (1.6)  on  fx,  the  map  fx  undergoes  a  generic  period  doulbling 
at  A  =  0.  In  fact,  there  is  a  neighborhood  I'  x  IF  C  x  R  of  (0,0)  G  IR  x  /R  such  that, 
for  each  value  of  A  G  I'F  for  which  X0  >  0,  there  is  a  unique  period  two  point  {o:x,P\)  of 
fx  with  o:x,Px  €  V.  We  say  that  the  bifurcation  is  supercntical  if  >  0  (the  fixed  point 
0  is  stable  for  /o)  and  subcritical  d  <  0  (the  fixed  point  0  is  unstable  for  /o)  (see  Figure 
1.5).  This  implies  that  the  period  two  point  is  Iwperbolic  stable  (resp.  unstable)  if  /3  >  0 
(resp.  /?  <  0): 


^>0 


f)<0 


Figure  1.5.  Bifurcation  diagram  for  fx 
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Theorem  1.1.  Suppose  that  S,  U  are  the  stable  and  unstable  regions  of  the  origin  as¬ 
sociated  with  the  Hopf  bifurcation  curve  and  that  fx  satisfies  (1-6).  Then  there  is  a 
neighborhood  V  of  (0,  0)  in  the  (A,  e)-plane  and  a  neighborhood  W  of  x  =  0  such  that,  if 
/3  >  0  (resp.  /?  <  Oj  and  (A,  e)  £  V,  then  there  is  a  periodic  solution  of  (1.4)  in  W  with 
period  2r(A,  c)  =  2  +  2c  +  0(|c|(|Al  +  |c|))  as  (A,  c)  (0,  0)  if  and  only  if  {\,  c)  €  F  n  t/ 
(resp.  (A,c)  G  Ffl  Sj.  Furthermore,  this  solution  is  unique  and 

(1)  »0ase— >0  uniformly  on  compact  subsets  of  lEl  \  JL  if  P  >  0, 

(2)  is  pulse-like  and  x^'^  —  p^—*0as€—*0  uniformly  on  compact  subsets  of 
JB.\^ifp<0, 

where  s^,  are  respectively  the  square  and  pulse  2-periodic  functions  defined  above.  If, 
in  addition,  fx(z)  =  -fx(-z),  then  x^‘^(i  +  t(A,  c))  = 

In  the  case  0  <  0,  the  meaning  of  pulse-like  is  that  the  solution  limits  to  constants 
on  the  integers.  However,  the  values  of  these  constants  exceed  the  values  of  the  correspond¬ 
ing  period  two  points  of  the  map.  Theorem  1.1  is  due  to  Chow,  Hale  and  Huang  (1992)  in 
the  supercritical  case  and  to  Hale  and  Huang  (1992a)  for  the  subcritical  case. 

It  is  possible  to  make  further  reasonable  conjectures  about  the  relationship  between 
the  map  fx  and  equation  (1-4)  for  e,  A  small.  More  specifically,  we  should  consider  the 
parameter  A  as  a  vector  and  obtain  a  generic  unfolding  of  a  codimension  g  singularity 
for  period  2  points  of  the  map  fx  and  then  obtain  a  relationship  between  the  bifurcation 
surfaces  in  A-space  for  period  two  points  and  the  bifurcation  surfaces  in  (e,  A)-space  for 
periodic  solutions  of  (1.4)  of  period  approximately  2. 

In  the  next  section,  we  outline  the  proof  of  Theorem  1.1.  In  Section  3,  we  given  some 
extensions  to  matrix  equations  and,  in  Section  4,  we  given  further  extensions  to  the  case 
where  we  have  differential  equations  coupled  with  difference  equations. 

There  are  results  concerning  the  existence  of  periodic  solutions  which  are  not  neces¬ 
sarily  small.  In  fx(x)  represents  a  negative  feedback,  xfx(x)  <;  0  for  i  51^  0  and  leaves 
an  interval  invariant  (Theorem  2.4,  Part  1),  there  exist  a  nontrivial  periodic  solution  for 
every  c  >  0,  A  >  0.  Under  some  additional  conditions,  the  limiting  profile  as  €  0  is 

a  square  wave  (see  Mallet-Paret  and  Nussbaum  (1986a)).  When  related  to  Theorem  1.1, 
these  conditions  correspond  to  the  case  0  >  0;  that  is,  the  orbit  is  stable.  No  global  results 
which  correspond  to  the  case  0  <  0  are  known. 

The  conjecture  above  relating  the  periodic  solutions  of  (1.4)  to  the  period  two  points 
of  the  map  fx  is  known  not  to  hold  in  general  if  the  solutions  are  not  required  to  remain  in 
a  small  neighborhood  of  the  origin.  Furthermore,  in  this  global  setting,  it  is  generally  not 
possible  to  associate  period  k  points  of  the  map  with  solutions  of  (1.4)  (see  Mallet-Paret 
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amd  Nussbaum  (1986b),  (1994)).  On  the  other  hand,  these  global  results  do  not  rule  out 
the  possibility  of  the  validity  of  the  conjecture. 

2.  Proof  of  Theorem  1.1.  We  now  give  an  outline  of  the  ideas  of  the  proof  of  Theorem 
1.1,  omitting  the  nontrivial  technical  details  although  they  have  independent  interest. 

.A.t  a  point  (Ao,co)  on  the  Hopf  bifurcation  curve,  there  are  two  purely  imaginary 
solutions  ±zu;o  of  (1.2)  and  the  remaining  solutions  have  negative  real  parts.  It  is  possible  to 
extend  the  classiceil  transformation  theory  (theory  of  normal  forms)  in  ODE  for  determining 
the  approximate  flow  on  the  center  manifold  of  (1.4)  at  (A,  e)  =  (Ao,€o)  corresponding  to 
these  purely  imaginary  eigenvalues.  If  we  assume  that  the  function  fx  is  given  by  (1.6) 
with  0^0  and  these  computations  are  performed,  then  the  stability  properties  of  the 
origin  under  the  mapping  fx  are  determined  by  0  and  the  origin  is  stable  (resp.  unstable) 
\i  0  >  Q  (resp.  0  <  0),  which  corresponds  to  the  generic  supercritical  (resp.  subcritical) 
bifurcation  of  the  map  /a  at  A  =  0.  Once  this  normal  form  hcis  been  obtained,  then  the 
fact  that  the  solutions  of  (1.3)  close  to  ±ia;o  cross  the  imaginary  axis  from  right  to  left  as 
€  decreases  implies  that  there  is  a  generic  supercritical  (resp.  subcritical)  Hopf  bifurcation 
with  respect  to  c  at  (Ao,eo).  Therefore,  in  a  neighborhood  of  any  point  (Ao,co)  on  the 
Hopf  bifmcation  cvurve,  there  is  a  unique  periodic  orbit  in  the  region  S  (resp.  U)  il  0  >  Q 
(resp.  0  <  Q)  which  is  stable  (resp.  unstable).  Of  course,  we  do  not  know  if  this  orbit 
exists  and  is  unique  in  all  of  S  (resp.  17);  of  course,  in  a  small  neighborhood  of  the  origin. 
It  will  be  one  of  tasks  to  show  that  this  is  the  case. 

Normal  form  theory  for  retarded  functional  differential  equations  is  an  interesting  md 
nontrivial  subject  in  itself  md  we  refer  the  reader  to  the  papers  of  Faria  and  Magalhaes 
(1991),  (1992)  for  the  theorj’^  as  well  as  applications. 

We  remark  that  we  could  have  obtained  the  same  information  about  the  generic  Hopf 
bifurcation  by  using  the  method  of  Liapunov-Schmidt  to  obtain  the  bifurcation  function 
for  the  periodic  orbits  in  a  neighborhood  of  the  origin.  For  a  discussion  of  this  approach, 
see,  for  example,  Stecli  (1979),  (1985). 

Returning  to  our  proof,  we  next  show  the  existence  of  the  periodic  orbit  in  5  (resp. 
U)  and  then  determine  the  limiting  profile  of  the  periodic  solution  as  c  — ♦  0.  Since  we 
expect  to  have  some  fast  changes  in  the  shape  of  the  periodic  solution  near  the  integers,  we 
enlarge  the  system  by  introducing  new  variables  which  represent  separately  the  solution 
on  the  interval  (0, 1  +  re)  and  [1  +  re,  2  +  2r£).  Our  tranformation  also  will  change  the 
bifurcation  problem  in  the  ‘bawl’  parameters  (A,  e)  (bad  because  we*  must  consider  c  >  0 
and  cannot  consider  it  in  a  neighborhood  of  zero)  into  a  bifurcation  problem  with  ‘good’ 
parameters  (A,  r)  (good  in  the  sense  that  (A,r)  can  be  considered  in  a  full  neighborhood 
of  some  point).  To  accomplish  this,  we  introduce  some  scalings  which  were  originally 
proposed  by  S.-N.  Chow  several  years  ago.  We  suppose  that  (1.4)  has  a  periodic  solution 
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z{t)  with  period  2  +  2rc  and  we  let 


(2.1) 


Wi{t)  =  z(-crt),  102(0  =  +  1  +  tr) . 


Since  2{t)  has  period  2  +  2re, 

(2.2) 


we  see  that 

W2{t)  =  z{-er(t  +  1)  -  1) 
W2{t  —  1)  =  z{—ert  —  1) . 


If  we  use  (2.1)  and  (2.2)  in  (1.4),  we  deduce  that 

2^  ^i(0  =  »*«^i(0  -  rfx{tV2{t  -  1)) 

W2{t)  =  rw2{t)  -  rfxiwi{t  -  1))  . 


This  equation  is  now  independent  of  £.  We  remark  that,  if  (u>i,u;2)  is  a  solution  of  (2.3), 
then  (102,101)  also  is  a  solution  of  (2.3). 

We  next  determine  the  approximate  value  of  the  constant  r  in  the  period  2  +  2re.  This 
is  obtained  by  considering  the  linear  variational  equation  around  the  zero  solution  of  (2.3) 
for  A  =  0, 

^2  4^  ^i{t)  =  rwi{t)  +  rw2{t-l) 

102(0  —  ^^2(0  +  rwi(t  —  1) . 


The  eigenvalues  of  (2.4)  are  the  roots  of  the  characteristic  equation, 
(2.5)  det(;iJ  —  rLe^'I)  =  ()ti  —  r)^  —  =  0 , 


where 


(2.6) 


L^p  —  v?(0)  + 


1 

0 


v?(-l). 


The  left  hand  side  of  Equation  (2.5)  always  has  /i  =  0  as  a  zero.  It  is  a  simple  zero  if  r  ^  1 
and  it  is  a  double  zero  if  r  =  1.  Bifurcation  from  a  simple  zero  can  never  lead  to  einy  periodic 
orbits.  Therefore,  we  are  forced  to  take  r  =  1  in  the  first  approximation.  For  r  =  1,  the 
remmning  eigenvalues  of  (2.4)  have  negative  real  parts.  If  we  let  r  =  1  +  to  =  (lui,  u;2), 
where  5  is  a  small  p£irameter,  then  (2.3)  can  be  written  as 


(2.7) 


w{i)  =  Lwt  +  SLwt  -  Fx,6(tvt) , 


where 

(2.8) 


Fx,s(<p)  =  (1  +  5) 


V>2(— 1)  +  /a(¥22(~1)) 

<p=l 

y>2 

1 
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and  Wi{6)  =  ‘U){t  +  6)  for  —1  <  0  <  0. 

We  now  consider  Equation  (2.7)  as  a  perturbation  of  the  linear  equation 

(2.9)  0(t)  =  Lvt . 

Of  course,  we  will  consider  (2.7)  with  initial  data  in  the  space  C  =  C([-l,  0],  E?).  Since 
the  cliaracteristic  equation  for  the  linear  part  of  (2.7)  for  (A,  ^)  =  (0,  0)  has  a  zero  as  a 
root  of  multiplicity  two,  we  know  that  the  small  periodic  orbits  of  (2.7)  will  lie  on  a  two 
dimensional  center  manifold  which  is  tangent  to  the  subspace  generated  by  generalized 
eigenvectors  ,  (^2  associated  with  the  eigenvalue  zero  of  (2.8). 

To  obtain  a  center  manifold,  we  use  the  usual  decomposition  theory  for  linear  systems 
(see  Hale  and  Verduyn-Lunel  (1993))  writing  C  =  P  ®  where  P  is  the  generalized 
eigenspace  of  0  and  Q  is  a  subspace  complementary  to  P  and  invariant  under  the  semigroup 
defined  by  (2,9),  If  we  let  Wt  =  +  Z2C2  +  Wt^  where  Wt  G  Q,  then  a  center  manifold  is 

given  by  a  function  ^(^1,^2)  which  vanishes  at  (0,0)  and  is  defined  by  an  integral  equation 
which  is  obtained  from  the  variation  of  constants  formula.  To  solve  our  problem,  we  need 
to  know  the  specific  form  of  the  vector  field  on  the  center  manifold.  This  is  a  nontrivial 
computation  very  similar  in  spirit  to  those  that  are  involved  in  nomal  form  theory.  If  we 
perform  these  computations,  we  see  that  the  approximate  flow  on  the  center  manifold  is 
given  by  the  system  of  ordinary  differential  equations 

(2.10) 

22  =  -2l 


up  through  terms  of  order  (S  +  A)^|2|  +  \z\^. 

Because  of  the  symmetry  in  (2.3),  it  can  be  shown  that  the  periodic  orbits  of  (2.10) 
which  encircle  the  origin  and  have  period  >  2  are  in  one-to-one  correspondence  with  the 
periodic  solutions  of  (1.4)  of  period  2  +  2r€.  Of  course,  in  this  statement,  we  assume  that 
c  is  small  as  well  as  the  periodic  orbits  that  are  being  considered.  With  this  observation, 
it  is  sufficient  to  look  for  these  types  of  periodic  solutions  of  (2.10)  in  a  neighborhood  of 
the  origin  regarding  (A,  5)  as  the  bifurcation  pzirameters.  We  remark  that  this  problem 
is  different  from  the  usual  bifurcation  problem  for  periodic  orbits  since  we  are  concerned 
only  with  those  periodic  orbits  that  encircle  the  origin. 

We  can  apply  the  theory  of  normal  forms  to  (2.10)  to  make  a  nonlinear  change  of 
variables  in  (2.10)  which  is  close  to  the  identity  and  arrive  at  an  equivalent  equation 

^2  ~~  ~  4)9ziZ2 

Z2  =  -Zi 
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up  through  terms  of  order  (S  +  A)^l2|  +  l^l^. 

The  complete  unfolding  of  the  singularity  in  (2.11)  is  known  (see  Takens  (1974a), 
(1974b),  Carr  (19S1)).  The  case  /3  >  0  has  a  simple  structure  for  the  periodic  orbits  (they 
all  encircle  the  origin)  whereas  the  case  <  0  is  very  complicated  (periodic  orbits  for  some 
values  of  the  parameters  encircle  points  other  than  the  origin).  As  remarked  earlier,  we 
are  interested  only  in  those  that  encircle  the  origin.  This  means  that  the  unfolding  theory 
does  not  help  in  this  problem.  However,  the  techniques  used  in  this  theory  can  be  adapted 
to  discuss  the  above  problem  and  to  complete  the  proof  of  Theorem  1.1.  As  we  will  see, 
the  computations  are  very  complicated  since  the  e.vpression  c,  which  determines  the  period 
of  the  solutions  of  (1.4),  involves  Abelian  integrals. 

Let  us  be  a  little  more  specific  about  how  this  is  done.  We  consider  only  the  case 
where  ,5  <  0  since  this  is  by  far  the  most  complicated.  To  analyze  the  periodic  solutions 
of  (2.11),  it  is  convenient  to  rescale  variables 


A  =  — h  =  /iS, 


ui(i)  = 


,  t 
\/2/r 


tioit)  = 


,  ‘  ' 


in  (2.11)  to  obtain  the  equivalent  system 


(2.12) 


where  7 


ii\  =  uo 

Uo  —  QUi  +  7^2  “  —  '2\/2flU2Ul  +  0{fx~[ui  4-  1x2]^) , 

—  \/2{6  —  jO'/x).  For  /X  =  0,  S  =  0,  we  obtain  the  conservative  system 


(2.13) 


with  the  first  integral 


lil  =  U2 

1x2  =  atxi  —  Uj 


(2.14) 


TT  /  \ 


Equation  (2.12)  has  the  same  form  as  Equation  (4.4.1)  in  Carr  (1981,  p.64)  and  so  we  can 
use  his  technique  with  the  first  integral  (2.14)  to  find  conditions  on  /x,S  to  ensure  that 
there  are  periodic  solutions. 


Figure  2,1.  Phase  portrait  for  the  conservative  system  for  o:  >  0. 
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The  phase  portrait  for  the  Hamiltonian  system  (2.13)  is  shown  in  Figure  2.1  for  a  >  0. 
There  are  periodic  orbits  encircling  the  equilibrium  points  (±>/a>0)  However,  these  are  of 
no  interest  to  us  since  they  do  not  encircle  the  origin.  There  are  two  homoclinic  orbits  to 
the  origin  which  form  a  figure  eight.  There  also  are  periodic  orbits  encircling  the  figure 
eight  and  these  are  of  interest  to  us. 

If  a  <  0  there  is  only  one  equilibrium  point  of  (2.13)  and  all  orbits  are  periodic  and 
encircle  the  origin.  We  are  interested  in  all  of  these. 

To  understand  the  effect  of  the  perturbation  terms  in  (2.12),  we  compute  the  derivative 
of  the  function  H(a,ui,U2)  along  the  solutions  of  (2.12)  to  obtain 

(2-15)  H(,a,uu<i2)=l4-  2 + . . . , 


where  . . .  represents  higher  order  terms.  We  want  to  find  a  periodic  solution  near  to  the 
periodic  solution  of  the  conservative  system  that  passes  through  the  point  (c,0),  c  >  0  or 
equivalently  through  the  point  Q  =  (0,(y  -  on  the  U2-axis  (see  Figure  2.1).  If 

—  (i‘i(^»Q)5W2(Q))  is  the  solution  of  (2.12)  through  the  point  Q,  then,  for  /i  and  S 
sufficiently  small,  there  are  constants  t*_  <t\  such  that  t/2(<l,Q)  =  0,U2(t^,Q)  =  0. 

It  is  now  possible  to  show  that  the  symmetry  in  (2.3)  implies  that  u{t,  Q)  is  periodic  if  and 
only  if  H(a,ui(tl,Q),0)  =  if(a, Ui(t!j.,  Q),  0).  If  v?e  perfonn  an  integration  in  (2.15),  we 
deduce  from  this  remcirk  that 

(2-16)  — \/2(6  - -o/i)  =  2\/fP(Q;,  c)/x  +  . . . , 

where  . . .  denotes  higher  order  terms  and 


(2.17) 


Jo(a,c)  =  jf  (aV-^  +  y-ac^/^du,, 


where  c  >  -y/la  if  a  >  0  and  c>  0  if  a  <  0. 
For  a  >  0,  c  >  V^,  let 


I(a,  c)  =  {  (a,  c) :  a  G  [-a,  a],  c  >  c>  ifa>0,  c>c>0ifa<0}. 


We  can  now  apply  the  Implicit  Function  Theorem  to  (2.16)  to  obtain  the  following 
result. 
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Lemma  2.1.  For  each  Rxed  a  >  0,  c  >  there  are  a  constant  p.  >  0  and  a  function 

r(a,  c,  fi)  =  2V^P(q-,  c)ii  + 

defined  for  0  <  ^  <  (a,c)  G  /(a,c),  such  that  (2.12)  has  a  unique  periodic  solution 

passing  through  Uj  =  0,  1^2  =  [cV2  -  qc2]^/2  if  and  only  if 

9 

(2.18)  -v/2(<5  -  ^ap)  =  T{a,  c,  p)  =  2v/2P(a,  c)p  +  0{p^) . 

Furthermore,  if  a  >  0,  then  the  periodic  solution  tends  to  the  pair  of  homocUnic  orbits 

(the  figure  eight)  of  (2.12)  as  c  — »  \/2a  and,  if  a  <  0,  the  periodic  solution  tends  to  zero 
as  c  — +  0. 

Formula  (2.18)  determines  <5  as  a  function  S^p)  of  p  so  that  there  is  a  periodic  solution 
of  (2.12)  with  initial  data  Q  =  Q(c)  =  (0,  [cV2  -  ac2]i/2).  Let  6  =  6(p)  and  let  u{t)  = 
(ui(t,  p,  a,  c),  U2{t,  p,  a,  c))  (resp.,  uo{t)  =  a,  c),  U2o(t,  a,  c)))  be  the  periodic 

solution  of  (2.12)  (resp.,  of  (2.13))  of  period  T(p,  a,  c)  (resp.,  To(a,  c))  passing  through 
Q  -  (0,  [c'‘/2  -  q'c2]i/2).  Since  u(t,  p,  a,  c)  -4  uo{t,  a,  c)  as  p 0,  we  have  T{p,  q,  c)  -4 
^0(0;,  c)  as  //  — »  0.  Furthermore, 

r(/^,  O!,  c)  =  To(a,  c)  +  0{p) . 

From  the  above  scaling,  this  periodic  solution  leads  to  a  periodic  solution  of  (2.11)  of 
period  T(p,  a,  c)/ \/2p,  which  in  turn  gives  a  periodic  solution  of  (2.4)  of  period  2  +  2re 
if  and  only  if  2  +  2r6  =  er{T(p,  a,  c)/^.p),  \  ^ -p-^  and  r  =  I  +  pS,  where  6  is  given  by 

(2.18) .  Therefore, 

(2.20)  c  =  e{p,  a,  c)  =  - _ 

[r(/x,  a,  c)  -  2y/2p][l  +  /x2(^  -  2P(a,  c))]  ‘ 

Here,  we  omit  the  terms  of  order  O(p^)  since  they  play  no  essential  role. 

Lemma  2.2.  Let  c  >  \/25  be  given,  let  p  be  defined  as  in  Lemma  2.1  and  let  e{p,  a,  c) 
be  given  by  (2.20).  Then,  for  each  Rxed  p  €  (0,  p]  and  a  €  [0,  a],  the  function  e(p,  a]  c) 
IS  monotone  mcreasmg  for  c  €  (\/2a,  c]  provided  that  €(p,  a,  c)  >  0.  In  addition, 

lim  e(p,a,c)  =  0. 

Wc  remark  that  the  periodic  orbit  {  u(f ,  p.  a,  c) ,  t  e  R  )  tends  to  a  pair  of  homoclinic 
orbits  (the  figure  eight  of  (2.12))  as  c  v^.  Therefore,  T(p,  a,  c)  oo  as  c  \/5^, 
which  implies  the  relation  (2.21). 
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The  proof  of  this  lemma  is  the  most  difficult  part  of  the  proof  because  of  the  difficulties 
in  understanding  how  ratios  of  the  Abelian  integrals  in  (2.17)  depend  upon  the  parameter 
c.  We  refer  to  Hale  and  Hueing  (1992a)  for  the  details. 

Lemma  2.2  completes  the  proof  of  Theorem  3.1. 


3.  Hybrid  systems.  In  applications,  we  often  encounter  equations  which  are  gener¬ 
alizations  of  (3.1)  and  consist  of  differential  difference  equations  coupled  with  difference 
equations.  We  refer  to  such  systems  as  hybrid  systems  and  the  form  that  we  consider  is 
the  following: 

ex(t)  =  Mxit),yit)) 

y(0  =  9\{x{t),y{t),x{t  -  l),y(<  -  1)) , 

where  e  >  0,  A  are  small  real  parameters,  i  €  R"*,  y  6  R”  are  vectors  and  the  functions 
fx{^-.y)  —  /(^iVi  A)  and  gx{x,y,z,w)  =  y(x,y,  2,«;,  A)  are  smooth  vector  valued  functions 
which  vanish  for  all  variables  equal  to  zero. 

This  system  includes  (1.4)  by  letting  the  function  gx  be  x{t  -  1)  and  the  function  fx 
be  replaced  by  -x{t)  +  fxiy{t)).  System  (3.1)  includes  also 


(3.2) 


€i(i)  +  x(t)  =  fx(x{t  -  l),x(t  -  2)) 


by  letting  y  be  a  two  vector  y,(t)  =  x(<  -  l),y2(t)  =  yi(t  -  1)  =  x{t  -  2).  In  this  way, 
we  are  able  to  consider  situations  in  which  the  map  defined  by  e  =  0  is  a  two  dimensional 
map;  for  example,  the  Henon  map. 

In  models  of  transmitted  light  through  ring  cavities  with  several  chcimbers,  Vallee, 
Dubois,  Cote  and  Delisle  (1987),  Vallee  and  Marriott  (1989),  have  proposed  the  following 
model: 

(3-3)  ^  ~  hx(z{t  -  1)) , 

where  each  >  0  is  a  small  parameter  and  2  is  a  scalzir. 

We  can  transform  this  to  an  equivalent  matrix  system 

^iii(i)  +  ^i(0  =  3^2(0 


(3.4) 


— l^m  — 1(0  d”  ®m— l(f)  — 

^m®m(0  d"  ^fn(f)  —  hx(^Xi(i  —  1))  . 


If  we  assume  that  Sj  as  Sj  =  eaj  j  =  1,  2, 
then  we  obtmn  the  matrix  equation 


m,  where  each  a,  >  0,  i  =  1, 2, . . . ,  m. 


(3.6) 


£i(t)  +  Ax(<)  =  A/A(i(t-l)), 
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where  A  and  fx  are  given  by 


(3.6) 


■^1 

-ai 

0 

0 

0  • 

‘hx(xi)- 

0 

02 

-02 

0 

0 

hxixi) 

• 

• 

• 

... 

• 

= 

• 

0 

0 

0 

...  Ofji^i 

-“m-l 

hx(xi) 

.  0 

0 

0 

0 

-hxixi). 

System  (3.5)  is  a  special  case  of  (3.1)  if  we  put  y{t)  —  x{t  —  1). 

It  is  of  interest  to  consider  (3.5)  for  a  general  function  fx  and  a  general  matrix  A, 
assuming  only  that  A~^  exists.  Such  a  system  is  a  special  case  of  the  system 


ex{t)  +  Ax{t)  =  Afxiyit)) 
y(0  =  5a(x(<  -  1),  y(t-l)), 


where  i  6  1R"\  y  G  R"  are  vectors,  the  m  x  m  matrix  A  has  zui  inverse  and  the  functions 
fx{y)  and  gx{^,  y)  are  smooth  vector  valued  functions.  Equation  (3.7)  has  been  used  by 
Ikeda  (1979),  Ikeda,  Daido  and  Akimoto  (19S0)  as  a  model  of  a  ring  cavity  containing  a 
nonlinear  dielectric  medium  for  which  a  part  of  the  transmitted  light  is  fed  baclv  into  the 
medium. 

The  equation  (3.1)  zirises  also  in  the  theory  of  transmission  lines.  If  the  lines  are 
lossless  and  described  by  the  telegraph  equations  with  the  boundary  conditions  for  the 
circuitry  between  the  lines  reflecting  Kirchoff’s  laws,  it  has  been  known  for  a  long  time 
that  the  flow  can  be  described  bj'  an  equivalent  set  of  neutral  delay  differential  equations 
(see,  for  example.  Hale  and  Verduyn-Lunel  (1993)  for  a  discussion  and  references).  Many  of 
these  same  problems  ^so  can  be  written  in  the  form  (3.1).  For  example,  by  an  appropriate 
change  of  coordinates,  the  equation  studied  by  Shimura  (1967)  for  a  transmission  line  with 
a  tunnel  diode  and  a  lumped  parallel  capacitance  can  be  rewritten  as 

,3  gs  ci(t)  =  y(t)  -  g(x(t)) 

y(t)  =  a  +  Ky{t  -  1)  —  x{t)  —  Lx{t  —  1)  , 

where  (x(t),  y(t))  represent  the  voltage  and  current  at  one  end  of  the  line,  all  constants  are 
positive  and  represent  physical  parameters.  Under  reasonable  assumptions  in  the  model, 
the  pzirameter  e  can  be  considered  to  be  very  small.  In  the  paper  of  Shimura  (1967),  several 
wave  forms  were  observed  numerically  which  compared  reasonably  well  with  experimental 
results.  Some  of  these  wave  forms  were  very  similar  to  square  waves. 

As  remarked  before,  in  many  situations,  for  fixed  €0  >  0,  there  is  a  A*(€o)  (which  we 
assume  to  be  positive  for  definiteness)  such  that  (3.1)  undergoes  a  generic  Hopf  bifurcation 
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at  (e,  A,x,y)  =  (co,  A*(co),0,0)  to  a  periodic  solution  Let  us  assume  that  there 

is  a  neighborhood  U  of  (e,A)  =  (0,0)  such  that  exists  for  all  (c  >  0,A)  in  U. 

Our  objective  remains  to  understand  the  behavior  of  the  profile  of  y*^)  as  c  -+  0.  We 
are  not  able  to  do  this  in  the  general  context  described,  but  we  can  say  something  if  we 
impose  more  conditions  on  the  functions  f\,gx. 

It  is  to  be  expected  that  the  limiting  profile  is  in  some  way  related  to  the  equation 
obtained  by  putting  c  =  0  in  (3.1).  For  e  =  0,  we  suppose  that  the  resiilting  equation  (3.1) 
defines  a  map  on  IR"'  x  IR" : 

(3.9)  (x,y)  6R"’ xR"H-»rA(x,y)elR'"  xlR" 

for  which  the  origin  is  stable  for  A  <  0  and  unstable  for  A  >  0.  Let  us  also  suppose  that 
Tx{x,y)  undergoes  a  generic  period  doubling  bifurcation  at  (x,y,  A)  =  (0,0,0)  with  the 
period  two  points  being  djx  €  IR"*  x  IR",  j  =  1,2.  If  the  bifurcation  is  supercritical,  we 
can  define  the  square  wave  2-periodic  function  and,  if  the  bifurcation  is  subcritical,  we  can 
define  the  pulse  wave  2-periodic  function  as  in  Section  1. 

Under  some  conditions  on  the  functions  fxiQXt  the  generic  period  doubling  bifurcation 
of  Tx{x,y)  leads  to  a  generic  Hopf  bifurcation  in  (3.1)  which  is  supercritical  (resp.  sub- 
critical)  if  the  period  doubling  bifurcation  is  supercritical  (resp.  subcritical).  The  natural 
question  as  before  is  the  following: 

Is  U  possible  that  the  limiting  profile  of  the  periodic  solution  (x*_)^,  y*;^)  obtained  through 
the  Hopf  bifurcation  is  either  the  square  wave  or  pulse  wave? 

In  a  later  section,  we  present  general  results  of  Hale  and  Huang  (1994)  for  which  this 
is  true  for  equation  (3.7).  Preliminary  computations  indicate  that  similar  results  also  hold 
for  systems  of  the  form  (3.8). 

4.  Hopf  bifurcation  in  a  special  hybrid  system.  In  this  section,  we  give  conditions 
under  which  there  will  be  a  generic  first  Hopf  bifurcation  for  the  hybrid  system 

ei(f)  +  Ax{t)  =  Afxiyit)) 

y(0=5A(x(<-l),  y(<-l)), 

where  x  €  IR"*,  y  €  IR"  are  vectors,  the  m  x  m  matrix  >1  has  an  inverse  and  the  functions 
fxiv)  and  gx{x,  y)  are  smooth  vector  valued  functions. 

We  malce  the  following  hypothesis: 

(Hi)  exists,  /;^(0)  =  0,  yA(0,0)  =  0. 

(Hi)  is  natural  since  we  want  to  consider  bifurcation  firom  the  origin  and  we  also  want  a 
map  to  be  defined  for  e  =  0. 
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If  we  introduce  the  notation, 


Aj(A)  =  C,A(0).  B,(A)  =  Z).ja(0,0),  Bj(A)  =  D,!,a(0,0) 

C(A)=B,(A),4,(A)  +  B,(A), 

then  the  linear  variational  equation  of  (4.1)  about  the  origin  is 


ex{i)  4-  Ax{t)  =  A42(A)y(t) 
y(0  =  Bi{X)x{t  -  1)  +  B2{X)y(t  -  1) , 


for  which  the  characteristic  matrix  is 


(4.4) 


A(A,€,/i) 


efilm+A  --4>l2(A) 

-Bi(A)e-^  In-B2{X)e-*^ 


We  let  (t{C)  denote  the  spectrum  of  a  square  matrix  C  eind  let  Bp  =  {r  €  C :  |z|  <  p}. 
Our  next  hj'pothesis  is 

(H2)  <^(^2(0))  cBp,  p<  1. 


Without  hypothesis  (H2),  it  is  not  too  difficult  to  show  that  there  would  be  solutions  of 
(4.4)  which  either  accumulate  or  lie  on  a  line  in  the  complex  plane  with  real  part  >  0.  If 
this  were  the  case,  we  would  not  be  able  to  find  a  simple  Hopf  bifurcation. 

It  is  convenient  to  introduce  the  following  definition.  We  say  that  a  curve  F  in  the 
(A,e)-plane,  e  >  0,  is  a  Hopf  Bifurcation  Curve  of  (4.3)  if  there  is  an  e*  >  0  and  a 
continuous  function  A  =  A(€),0  <  e  <  e*,  A(e)  0  as  e  — +  0  such  that,  if  Ff  =  {(A,  e)  : 
A  =  A(e),  e  G  (0,  e*)},  then,  for  any  (Aq,  co)  ^  F<- ,  there  are  two  purely  imaginary  solutions 
±i73o  of  the  characteristic  equation  det  A(Ao,€o,p)  =  0  and  the  remaining  solutions  p 
satisfj'  Re  /X  7^  0.  We  say  that  a  Hopf  Bifurcation  Curve  is  a  Generic  Hopf  Bifurcation 
Curve  wWi  respect  to  e  if,  for  fixed  Aq,  the  two  eigenvalues  /x(Ao,  e),  /2(Ao,  e),  /x(Ao,  Co)  =  t^Oj 
satisfy  Re5/x(Ao,eo)/5e  <  0.  This  type  of  transversal  crossing  of  the  imaginary  axis  of 
the  eigenvalue  p(Xo,e)  implies  that  there  will  be  a  Hopf  bifurcation  with  respect  to  c  at 
Cq.  We  say  that  a  Hopf  Bifurcation  Curve  is  the  First  Hopf  Bifurcation  Curve  if,  for 
each  (Ao,€o)  €  Ffj,  all  eigenvalues  p  corresponding  to  the  parameters  e  >  cq?  A  =  Aq, 
have  Re  <  0.  The  Generic  First  Hopf  Bifurcation  Curve  with  respect  to  e  is  a  Generic 
Hopf  Bifurcation  Curve  and  also  a  First  Hopf  Bifurcation  Curve.  The  Generic  First  Hopf 
Bifurcation  Curve  with  respect  to  e  is  the  most  interesting  because  there  is  a  transfer  of 
stability  of  the  origin  at  e  =  £©;  that  is,  the  origin  is  stable  for  e  >  €o  and  unstable  for 
e  <  Co-  From  the  physical  origins  of  the  problem,  this  is  natural  because  we  expect  that  the 
origin  is  stable  for  large  c  (by  a  change  of  time  scale,  this  is  small  delay)  and  to  eventually 
become  unstable  for  small  c  (large  delay).  The  Generic  First  Hopf  Bifurcation  Curve  with 
respect  to  e  represents  the  first  change  in  the  stability  properties  of  the  origin. 
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Huang  and  Hale  (1994a)  have  pven  conditions  on  the  coefficients  in  (4.2)  for  the 
existence  of  the  Generic  First  Hopf  Bifurcation  Curve. 

In  addition  to  (H2),  assume  that 


(H3) 


—(1  +  A)  €  <t(C(A))  is  a  simple  eigenvalue, 

\  {~(1  +  ‘^)}  C  Bp  with  p  <  1. 

If  (H3)  is  satisfied,  then  we  can  introduce  a  change  of  variables  in  y  to  obtain 


(4.5)  C(A)  = 

We  also  introduce  the  notation 

(4.6)  -Bi{0)A-^A2i0)  =  S  = 


-(1  +  A)  0 

0  HoiX) 


,  a{Ho{X)  cBp,  p<  1. 


Sii  Si2 
S21  S22 


,  -Bi{0)A-^A2{0)  =  w  = 


Wn  Wi2 
W21  W22 


The  remaining  hypotheses  that  we  impose  are 


(H4)  -Ro  =  — •S'lj  +  2[Wii  +  5i2(/ +  ifo(O))  7^  0,  5ii  0. 


(H5) 


min{Re2: :  z  €  o{A)]  >  0, 


(H6) 


Det 


iOIm  +  A 
-51(0)6-- 


-A^2(0) 

In  -  52(0)e-- 


^  0,  for  0  >  0,0  <  V  <  27r, (^,t;)  ^  (0,7r) 


(H7) 


Det 


p/m  +  A 

.  -5i(0) 


-.4^2(0) 

In  -  52(0) 


7^  0,  for  all  p  G  C,  Re  p  >  0. 


One  of  the  main  results  of  this  section  is 

Theorem  4.1.  Under  the  assumptions  (Hl)-(H5),  the  hypotheses  (H6)  and  (H7)  are 
necessary  and.safEcient  for  the  existence  of  the  Generic  First  Hopf  Bifurcation  Curve  with 
respect  to  e.  Furthermore,  if  these  hypotheses  are  satisfied,  then  there  is  an  e*  >  0  such 
that  this  curve  is  given  as  =  {((A(€),  c),  e  €  (0,  €*},  where  A(e)  is  a  C^-function  having 
the  property  that  A(c)  0  and  is  given  approximately  by 

(4.7)  A(e)  =  +  o(€^) 


as  e  —*  zero. 


Formula  (4.7)  implies  that  the  Generic  First  Hopf  Bifurcation  Curve  with  respect  to 
€  in  the  (A,  €)-plane  is  the  graph  over  the  poative  (resp.  negative)  A-axis  if  Rq  >  0  (resp. 

Ro  <  0). 
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We  do  not  prove  this  theorem  and  refer  the  reader  to  Hale  and  Huang  (1994b). 

Let  us  see  if  the  hypotheses  of  Theorem  4.1  are  satisfied  for  some  of  our  examples. 
We  have  remarked  before  that  (1.4)  is  a  special  case  of  (4.1)  when  it  is  written  as 

(4  8)  «(0  +  x{<)  =  /A(y(i)) 

y(()  =  i((  -  1). 

Let  us  assume  that 

(4-9)  fx(x)  = —{1  +  X)x  +  ax^  +  bx^, 

where  0  =  ar  +  b  ^  0.  The  hneeir  variational  equation  about  a;  =  0  for  (4.8)  is  a  special  case 
of  (4.3)  with  A  =  1,  j42(A)  =  —(1  +  A),  Bi(A)  =  1,  B^iX)  =  0.  It  is  now  obvious  that  the 
hypotheses  (Hl)-(H5)  are  satisfied.  A  simple  computation  shows  that  hypotheses  (H6) 
is  equivalent  to  z(0  -  sin  v)  +  1  +  cos  u  5^  0  for  0  >  0,  0  <  u  <  2?:,  u  7^  tt.  This  is  clearly 
satisfied.  Also,  (H7)  is  equivalent  to  ^  +  2  7^  0  for  Re  /Li  >  0,  whicli  is  true.  Therefore, 
there  is  a  Generic  First  Hopf  Bifurcation  Curve  with  respect  to  e. 

Let  us  next  discuss  equation  (3.5),  (3.6).  We  show  that  all  of  the  hypotheses  are 
satisfied  by  analyzing  directly  the  characteristic  equation.  The  proof  is  not  difficult.  The 
characteristic  equation  for  the  linearization  about  a:  =  0  is  given  by 

(4-10)  E{e,  A,  fi)  =  {eamfJ-  +  1)  •  •  •  (earj/i  +  1)  +  (1  +  A)e~''  =  0. 

In  this  particular  situation,  the  determinant  in  (H6)  is  given  by 

Ei{d,v)  =  E{-,0,iv)  =  JJ(t6'a_,-  +  1)  +  e"’". 

j=i 

If  ^  >  0,0  <  V  <  2?:,  axe  such  that  £^1(^,1;)  =  0,  then  we  have  [  +  1)| 

is,  +  1)  =  1.  Since  aj  >  0,  we  must  have  6  =  0  and  hence  1  +  e"'” 

is,  u  =  TT.  Therefore,  the  condition  (H6)  is  satisfied. 

Next,  the  fimction  in  (H7)  has  the  form 

m 

c'2(/i) = + 1) + 1- 

i=i 

K  there  is  a  /i  =  u  +  iVyU  >  0,  such  that  C'2(u  +  tv)  =  0,  then  it  follows  that  1  = 
"b  1)  "b  ^  )•  Therefore,  we  must  have  u  =  v  =  0,  which  leads  to  the  assertion 
that  0  =  C2(0)  =  1  +  1,  which  is  a  contradiction.  As  a  consequence,  C2(/x)  7^  0  for  all 
€  C  with  Re/i  >  0;  that  is,  (H7)  is  satisfied. 


=  1;  that 
=  0;  that 
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We  can  now  show  that,  for  fixed  cq  >  0,  there  is  a  unique  A*(£o)  >  0  such  that  (4.10) 
heis  exactly  two  purely  imaginary  roots  and  the  remaining  ones  have  negative  real  parts  for 
(A,c)  =  (A(eo)>co)-  for  c  >  Co>  the  origin  is  asymptotically  stable  and,  for  0  <  e  <  coj  the 
origin  is  unstable.  In  this  way,  we  obtain  the  existence  of  a  Generic  First  Hopf  Bifurcation 
Curve  with  respect  to  c. 

After  we  know  that  a  Generic  First  Hopf  Bifurcation  Curve  with  respect  to  c  exists, 
the  next  step  is  to  determine  the  direction  of  bifurcation  for  the  nonlinear  equation.  To 
describe  the  result,  we  need  some  additional  notation.  It  turns  out  that  the  direction  of 
bifurcation  is  related  to  the  map 

(4.12)  y  6  R"  Tx{y)  =  gxUxivU)  €  R" 

obtained  from  (4.1)  for  e  =  0.  If  we  assume  that  C(A)  is  given  as  in  (4.5)  and  let  y  =  col 
(yiiy2)  €  R^  X  R"  =  col  {J^\Xi^2x)  €  R*  X  R”~^,  then  we  can  write 

(4  11)  =  -(1  +  +  yA(X)v2  +  +  Odlyif  +  !,J||!,2||  +  y\) 

+yfir.(A)  +  +  Odtef  +  vJIIkII  +  ||yf) 

35  (y,A)  -» (0,0). 

We  assume  that 

(H8)  R,  =  k2(0)[/-  fro(0)]-'ifi(0)  +  kliO)  +  *3(0)  #  0. 

The  following  result  is  due  to  Hale  and  Huang  (1994b). 

Theorem  4.2.  Suppose  that  (Hl)-(H8)  are  satisfied  and  r<-  is  the  Generic  First  Hopf 
Bifurcation  Curve  with  respect  to  e  given  in  Theorem  4.1.  If,  in  addition,  Rq  >  0,  then, 
for  each  fixed  Aq  =  A(eo),  cq  sufficiently  small,  system  (4.4)  has  a  generic  Hopf  bifurcation 
from  the  origin  with  respect  to  e  at  e*  which  is  supercritical  (resp.  subcritical)  if  Ri  >  0 
(resp.  Ri  <  0). 

The  proof  of  this  is  very  long  and  uses  results  from  a  very  general  Hopf  bifurcation 
theorem  of  Hale  and  Huang  (1994a).  The  idea  for  the  proof  is  simple  and  standard.  It 
is  first  observed  that  Rq  >  0  implies  that  the  origin  is  stable  for  e  >  c*  and  iinstable  for 
e  <  £*.  From  the  physical  point  of  view,  this  is  to  be  expected  and  therefore  the  assumption 
that  i2o  >  0  is  natural.  Also,  there  is  a  bifurcation  function  G(a,  e)  for  periodic  solutions 
near  the  origin,  where  a  is  approximately  the  amplitude  of  the  periodic  solution.  There  is  a 
periodic  solution  of  (4.1)  near  the  origin  if  and  only  if  there  are  (a,  c)  such  that  G(a,  e)  =  0. 
This  function  is  odd  in  a  and  it  is  shown  that  the  linear  term  in  a  has  the  same  sign  as 
—Rq.  Furthermore,  the  cubic  term  in  a  hzis  the  S2Lme  sign  as  eRi.  Therefore,  there  is  at 
most  one  periodic  solution  and  it  satisfies  the  properties  stated  in  the  theorem. 

The  following  result  is  in  Hale  and  Huang  (1994b)  and  is  a  standard  application  of 
the  method  of  Liapimov-Schmidt. 
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Theorem  4.3.  If  (Hlj,  (H3j  and  (HS^  are  satisfied,  then  there  is  a  generic  period  two 
doubling  bifurcation  of  fFx  at  (y,  A)  =  (0, 0).  More  precisely,  if  Ri\  >0,  then  there  are 
period  two  points  dix,d2x  of  Tx  such  that  ^x(dix)  =  d2x,^x(d2x)  =  dix,  which  is  stable 
for  i?i  >  0  (supercritical  bifurcation)  and  unstable  for  Ri  <  0  (subcritical  bifurcation). 

This  is  a  very  interesting  result  because  it  says  that,  if  we  have  iZo  >  0  and  the 
conditions  in  Theorem  4.2  for  Generic  First  Hopf  Bifurcation  with  respect  to  c  to  occur, 
then  also  there  is  generic  period  doubling  of  the  map  with  respect  to  A  and  the  direction 
of  bifurcation  of  the  Hopf  bifurcation  with  respect  to  e  is  the  same  as  the  direction  of 
bifurcation  of  the  period  doubling  bifurcation  of  the  map  with  respect  to  A. 

We  remark  that  the  converse  of  this  statement  may  not  be  true.  More  precisely, 
it  is  possible  for  the  map  in  (4.12)  to  generically  period  double  at  A  =  0  and  yet  the 
lineeirization  about  the  origin  in  the  differential  euation  may  not  possess  a  Generic  First 
Hopf  Bifurcation  Curve  with  respect  to  c. 


5.  Square  and  pulse  waves  for  hybrid  systems.  It  remains  to  show  that  the  periodic 
orbit  obtained  from  the  Hopf  bifurcation  can  be  extended  away  form  the  Hopf  bifurcation 
cur\'e,  is  unique  and  has  limiting  values  related  to  either  the  square  or  pulse  periodic 
functions  related  to  the  period  two  points  of  the  map.  We  want  to  apply  a  method  similar 
to  the  one  used  in  the  proof  of  Theorem  1.1;  namely  scaling,  center  manifolds,  etc. 

As  in  Section  2,  we  seek  periodic  solutons  of  (4.1)  with  a  period  2  +  2(ro  +  h)e,  where 
To  is  a  fixed  parameter  (to  be  determined  later)  which  depends  only  upon  the  matrices 
A, .42(0),  Bi(0),  ^2(0)  and  h  will  be  determined  as  a  function  of  e.  As  in  Section  2,  if 
is  such  a  solution  of  (4.1),  we  introduce  the  transformation 

(5  1)  =  a-'(-e(ro  +  fi)0,  "2(0  =  2:(-e(ro  +  h)<  +  1  + e(ro  +  h)) 

(0  =  y(-e(»'o  +  h)t),  V2(t)  =  y(-e(ro  +  h)t  +  1  +  e(ro  +  h)). 

Since  x(i)  and  y(t)  have  period  2  +  2(ro  +  /i)e,  we  see  that 


u(t  -  1)  =  x(-e(ro  +  h)t  -  1) 
v(t  -  1)  =  y(-e(ro  +  h)t  -  1)  . 


If  we  use  (5.1)  and  (5.2)  in  (4.1),  we  deduce  that 


“i(0  =  (**0  +  fi)Aui(<)  -  (ro  +  h)Afx(vi(t)) 
«2(0  =  (ro  +  h)Au2(t)  -  (ro  +  h)Afx(v2(t))  ' 
Vi(t)  =  gx(u2(i-l),V2(i-l)) 

V2(i)  =  9x(ui(t  -  1),  t;j(t  -  1)) . 
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We  consider  /i,  A  as  the  new  bifurcation  patrameters.  We  must  determine  the  constant 
ro  and  the  procedure,  as  for  the  simpler  case  in  Section  2,  will  be  to  insist  that  the 
linear  variational  equation  of  (5.3)  about  zero  has  a  double  eigenvalue  zero  with  no  other 
eigenvalues  on  the  imaginary  axis. 

If  we  use  the  notation  of  Section  4  and  let  u  =  col  (ui,U2),  v  =  col  (^1,^2), 


(5.4) 


A  = 


'A  o' 

A-,  — 

■^2(0) 

0 

0  .4 

,  — 

0 

^2(0). 

0  Bj 

Bj{0)  0  J  ’ 


i  =  1,2, 


then  the  linear  variational  equation  of  (5.3)  for  (h,  A)  =  (0, 0)  is  given  by 


(5.5). 


u{t)  =  roAu{t)  —  roAA2v{t) 
v(t)  =  Biu{t  —  1)  +  B2v{i  —  1). 


The  eigenvalues  of  (5.5)  are  the  solutions  of  the  characteristic  equation 


(5.6) 


det  A(/x,ro)  =  0,  A(^,ro) 


fjL-roA  roAA2 
-Bi€->^  I-B2e->‘ 


Because  of  (2.3),  (2.4)  and  the  symmetry  in  (5.5),  zero  always  is  an  eigenvalue.  We  impose 
conditions  on  the  coefficients  to  ensure  that  0  is  an  eigenvalue  of  multiplicity  two  and  there 
are  no  other  eigenvalues  on  the  imaginary  axis.  To  do  this,  we  will  need  hypotheses  (Hi), 
(H2),  (H4)  of  Section  4  as  well  as  the  following  one: 


(H9) 


det 


Si  1-4.42 
hn  -  526-' 


^  0  for  a;  G  IR  \  {0}, 


where  Sn  is  given  in  (4.6). 

The  follwoing  results  are  due  to  Hale  and  Huang  (1994b). 

Lemma  5.1.  If  (Hi),  (H2),  (H4)  and  (H9)  are  satisfied  and  tq  =  Sn,  then  ft  =  0  is 
an  eigenvalue  of  (5.5)  of  multiplicity  two  and  there  is  a  S  >  0  such  that  the  remaining 
eigenvalues  satisfy  |Re/x|  >  5  >  0  and  there  are  only  a  finite  number  of  eigenvalues  with 
positive  real  parts. 


Theorem  5.2.  If  (H1)~(H4),  (H.8)and  (H.9)  are  satisfied,  then  there  is  a  a  neighborhood 
U  of  (0,  0)  in  the  (A,  e)  plane  and  a  sectorial  region  S  in  U  such  that,  if  (A,  c)  €  U,  then 
there  is  a  periodic  solution  XX, (  of  (2.4)  with  period  2r(A,  «)  =  2  +  2Sn€  +  0(|e|(|A|  +  |€|)) 
^  ^)  (0,  0)  if  and  only  if  (A,  e)  €  S.  Furthermore,  this  solution  is  unique. 

We  do  not  give  a  proof  of  these  results,  but  simply  refer  the  readed  to  Hale  and  Huang 
(1994b).  We  only  make  a  few  remarks  explaining  in  some  detail  other  properties  of  the 
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solutions  in  the  sector  S.  Of  course,  the  sector  5  must  belong  to  the  set  e  >  0  in  the  (A,  e) 
plane.  It  actually  is  shown  that,  if  >0  (the  supercritical  case  of  period  doubling  of  the 
map)  and  i2o  >  0,  then  the  sector  S  C  { (A,  c) :  e  >  0,  A  >  0  }  and,  for  A  =  Ao  >  0,  fixed, 
the  set  {  c  ;  (c,  Ao)  €  S' }  is  an  interval  (0,  eo(Ao)),’ where 

Ro  ^  ' 

as  Ao 0.  For  any  c  G  (0,eo(Ao)),  the  periodic  solution  xao.c(0  approaches  a  square  wave 
2-periodic  function  as  e  — »  0. 

If  <  0  (the  subcntical  case  of  period  doubling  of  the  map),  the  sector  S  is  com¬ 
pletely  different  (as  was  the  case  in  Theorem  1.1)  and  the  periodic  orbits  have  a  different 
structure  as  e  -»  0.  The  sector  5  contains  points  (e.  A)  with  A  both  negative  and  positive 
More  precisely,  for  Ao  >  0,  the  set  {e  :  (e,  Ao)  €  5}  is  an  interval  (eo(Ao),  ^o(Ao)).  For 
Ao  <  0,  the  set  {e  :  (e,  Ao)  €  S)  is  an  interval  (0,  ao(Ao)).  For  any  e  G  (0,Qo(Ao)), 
the  unique  periodic  solution  i;,,,,(t)  becomes  pulse  like  as  €  0  in  the  foUowing  sense: 

the  periodic  solution  xx^A^)  has  the  property  that  0  as  e  0  uniformly  on 

compact  sets  of  (0,  1)  U  (1,  2).  In  the  pulse  like  solution,  the  pulses  in  the  solution  occur 
near  the  integers  and  are  opposite  in  direction.  However,  the  magnitude  of  the  pulse  near 
the  integers  exceeds  the  magnitude  of  the  corresponding  period  two  point  of  the  map. 

Notice  that  Theorem  5.2  and  Theorem  4.2  use  different  hypotheses.  We  can  obtain 
genenc  first  Hopf  bifurcation  with  respect  to  e  (which  imply  the  generic  period  doubling 
of  the  map)  with  different  assumptions  than  for  the  existence  of  the  square  and  pulse  like 
solutions  related  to  the  period  2  points  of  the  map.  The  differences  arise  in  the  discussion 
of  the  eigenvalues  of  the  linear  variational  equation  for  (4.1)  and  the  one  for  the  scaled 
equations  (5.3).  It  is  possible  to  give  an  example  for  which  the  zero  solution  of  (4.3)  is 
stable  and  the  zero  solution  of  (5.5)  is  unstable.  Naturally,  this  leaves  open  some  questions 
as  to  whether  we  have  obtained  the  best  possible  results. 


References 

Berre,  MX.,  R«sayre,E.,  Tallel,  A.  and  H.M.  Gibbs  (1986),  High-dimension  chaotic 
attractors  of  a  nonlinear  ring  cavity,  Phys.  Rev.  Lett  56,  274-277. 

(1981),  Applicaiions  of  Centre  Manifold  Theory,  Appl.  Math.  Sci  35 
Springer- Verlag.  ‘  ’ 

Chow  S.-N.,  Hale,  J.K.  ^d  W.  Huang  (1992),  From  sine  waves  to  square  waves  in 
delay  equations.  Proc,  Royal  Soc.  Edinburgh  120A,  223-229. 

Chow,  S.-N.  ^d  W  Hu^g  (1992),  Singular  perturbation  for  delay  differential  equa¬ 
tions  and  apphcations.  J.  Differential  Equations.  To  appear.  ^ 

^7-  (^^53),  Singularly  perturbed  delay  differential  equa¬ 

tions.  In  Coupled  Nonlinear  Oscillators  (Eds.  J.  Chandra  and  A.  Scott),  North  Holland. 


34 


Faria,  T.  and  L.  T.  Magalhaes  (1991),  Normal  forms  for  retarded  functional  differential 
equations  and  applications  to  Bogdanov  singularity.  J.  Differential  Equations.  To  appear. 

Faria,  T.  and  L.  T.  Magalhaes  (1992a),  Normal  forms  for  retarded  functional  dif¬ 
ferential  equations  with  parauneters  and  applications  to  Hopf  bifurcation.  J.  Differential 
Equations.  To  appear. 

Faria,  T.  and  L.  T.  Magalhaes  (1992b),  Restrictions  on  the  possible  flows  of  scalar 
retarded  functional  differential  equations  in  neighborhoods  of  singularities.  J.  Differential 
Equations.  Preprint. 

Gibbs,  H.M.,  Hopf,  F.A.,  Kaplan,  D.D.L.  and  R.  L.  Shoemaker  (1981),  Observation 
of  chaos  in  optical  bistability,  Phys.  Rev.  Lett.  46,  474-477. 

Glass,  L.  and  M.  Mackey  (1977),  Oscillation  and  chaos  in  physiological  control  systems, 
Science  191,  287-289. 

Hale,  J.K.  and  W.  Huang  (1992a),  Period  doubling  in  singulcirly  perturbed  delay 
equations.  J.  Differential  Equations.  To  appear. 

Hale,  J.K.  and  W.  Huang  (1992b),  Square  and  pulse  waves  in  matrix  delay  differential 
equations.  Dyn.  Sys.  Appl.  1,  51-70. 

Hale,  J.K.  and  W.  Huang  (1993),  Variation  of  constants  for  hybrid  systems  of  FDE. 
Proc.  Royal  Soc.  Edinburgh.  To  appear. 

Hale,  J.K.  and  W.  Huang  (1994a),  Hopf  bifurcation  analysis  for  hybrid  systems.  Proc. 
Ini.  Conf.  Diff.  Equations  in  memory  of  Siavros  Busenburg.  To  appear. 

Hale,  J.K.  and  W.  Huang  (1994b),  Periodic  solutions  of  singularly  perturbed  delay 
equations.  Preprint. 

Hade,  J.K.,  Infante,  E.F.  and  F.-S.P.  Tsen  (1985),  Stability  in  linear  delay  equations. 
J.  Math.  Anal.  Appl.  105,  533-555. 

Hale,  J.K.  amd  S.M.  Verduyn-Lunel  (1993),  Introduction  to  Functional  Differential 
Equations,  Springer- Verlag,  New  York. 

Ikeda,  K.  (1979),  Multiple- valued  stationary  state  and  its  instability  of  the  transmitted 
light  by  a  ring  cavity  system.  Opt.  Commun.  30,  357. 

Ikeda,  K.,  Daiido,  H.,  amd  O.  .A.kimoto  (1980),  Optical  turbulence:  Chaotic  behavior 
of  transmitted  light  from  a  ring  cavity,  Phys.  Rev.  Lett.  45,  709-712. 

Mackey,  M.  and  M.  Milton  (1990),  Comments  on  Theoretical  Biology  1,  299-327. 

Mallet-Paret,  J.  and  Nussbaum  (1986a),  Global  continuation  and  asymptotic  behavior 
for  periodic  solutions  of  a  differential  delay  equation,  Annali  Mat.  Pura  Appl.  145,  33-128. 

Mallet-Paret,  J.  and  Nussbaum  (1986b),  Global  continuation  and  complicated  trajec¬ 
tories  for  periodic  solutions  for  a  differential  delay  equation,  Proc.  Symp.  Pure  Math.  45, 
Part  2,  155-167,  Am.  Math.  Soc. 

Mallet-Paret,  J.  and  Nussbaum  (1993),  Multiple  transition  layers  in  a  singularly  per¬ 
turbed  differential-delay  equation,  Proc.  Roy.  Soc.  Edinburgh  123A,  1119-1134. 

Shimura,  M.  (1967),  Analysis  of  some  nonlinear  phenomena  in  a  transmission  line. 
IEEE  Trans.  Circuit  Theory  14,  60-69.  , 

Stech,  H.W.  (1979),  The  Hopf  bifurcation:  stability  result  and  application.  J.  Math. 
Anal.  Appl.  71,  525-546. 

Stech,  H.W.  (1985),  Nongeneric  Hopf  bifurcation  in  functional  differential  equations. 
SIAM  J.  Appl.  Math.  16,  1134-1151. 


35 


(1974a),  Forced  oscillations  and  bifurcations.  Comm.  S,  Math.  Inst. 
Rijksuniversiteit,  Utrecht. 

Takens,  F.  (1974b),  Singularities  of  vector  fields.  Puhl.  Inst.  Hautes  Etudes  Sci  43 
47-100. 

Vallee,  R.,  Dubois,  P.,  Cote,  M.  and  C.  Delisle  (1987),  Phys.  Rev.  A36,  1327. 

Vallee,  R.  and  C.  Marriott  (1989),  Analysis  of  an  N‘*‘-order  nonlinear  differential  delay 
equation,  Phys.  Rev.  A39,  197-205. 

Wazewska-Czyewska,  M.  and  A.  Lasota  (1976),  Mate.  Stosovvana  6  (1976),  25-40. 


36 


Part  3,  Small  Delays  Can  Make  a  Difference 


1.  Difference  Equations.  In  Part  2,  we  have  discussed  in  some  detail  the  effects  of 
large  delays  on  the  dynamics  of  retarded  delay  differential  equations.  For  the  problems 
considered,  we  have  seen  that,  under  reasonable  assumptions,  the  limiting  dynamics  as 
the  delay  approached  infinity  presented  no  surprises,  at  least  locally  near  an  equilibrium 
point.  The  limiting  dynamics  was  determined  by  the  local  dynamics  of  a  map. 

If  we  consider  a  retarded  delay  differential  equation  (or  even  more  generally  a  func¬ 
tional  differential  equation)  for  which  the  delays  are  small,  then  it  is  possible  to  prove 
that  the  limiting  dynamics  is  determined  by  the  ordinary  differential  equation  obtained 
by  putting  all  of  the  delays  equal  to  zero  (see  Kurzweil  (1970),  (1971)  and  a  more  com¬ 
plete  discussion  in  Hale,  Magalhaes  and  Oliva  (1984)).  In  such  a  situation,  it  is  fair  to 
say  that  small  delays  are  unimportant.  For  retarded  functional  differential  equations,  the 
same  general  remarks  are  valid  if  we  make  small  changes  in  the  delays  with  the  change  not 
necessarily  occuring  around  zero. 

In  this  section,  we  devote  our  attention  to  similar  problems  for  neutral  differential 
difference  equations.  These  are  equations  for  which  the  derivative  of  the  solution  also 
occurs  with  a  delay.  Special  cases  of  such  equations  of  course  would  by  difference  equatons 
since  we  can  differentiate  to  obtain  a  neutral  differential  difference  equation.  In  the  next 
section,  we  will  see  also  that  these  problems  are  very  closely  related  to  the  control  of 
hyperbolic  PDE  when  the  control  function  is'  applied  on  the  boundary  with  a  time  delay. 

For  any  operator  T  on  a  Banach  space,  we  let  a(T)  denote  the  spectrum  of  T,  r(a(<)) 
denote  the  radius  of  the  spectrum  and  r^(a{T))  denote  the  radius  of  the  essential  spectrum. 
We  begin  with  an  example  of  a  diffei'ence  equation. 

Example  1.1  {Uniform  asymptotic  stability  may  be  destroyed  by  small  changes  in  the 
delay).  The  observations  in  this  example  are  essentially  due  to  Melvin  (1974).  Suppose 
that  Co  >  0,  0  <  ri  <  rj  <  Co,  are  constants  and  consider  the  difference  equation 

3:(t) -}- ^a:(t  -  ri ) -f  ix(f  -  ra)  =  0. 

Let  C  =  C([— co,0];R)  and  let 

C'o  =  C'o,(ri,ri)  =  {v’  €  C  ;  v>(0)  -1-  -(p{-ri)  +  i^(-r2)  =  0). 

For  any  V?  6  Co,  there  is  a  unique  solution  x(t,yi)  of  (1.1)  which  is  defined  for  all  t  G  R. 

If  we  define  (5'(ri,rj)(0v>](^)  =  x{t  -1-  6,(p),  6  €  I-Co,0],  then  S(ri,r{){t),  <  >  0,  is  a  C°- 
semigroup  on  Co,(r,,r,)- 

We  need  the  foUowing  well  known  result  (see,  for  example,  Hale  and  Verduyn-Lunel 

(1993)).  The  proof  is  nontrivial  and  makes  extensive  rise  of  properties  of  the  Laplace 
trarrsform. 
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Theorem  1.1.  »'(<^(%,,r,)(l))  =  e“,  where  a  is  defined  by 
(1-2)  a  =  a(rj,r2)  =  sup{ReA  :  Ao,(r,,rj)(A)  =  0}, 

We  now  compute  r(<^(-S’(r,,rj)(l))  for  some  particular  values  of  (ri,r2).  For  (ri,r2)  = 
(1,2),  Ao,(i,2)('^)  is  a  quadratic  function  in  e~^  and  it  is  ezisy  to  show  that  the  solutions 
of  (1.2)  satisfy  =  1/2  and  so  Re  A  =  -(ln2)/2.  Thus,  r(a(S(i,2)(l))  =  l/y/2  <  1, 

which  implies  that  the  zero  solution  of  (1.1)  is  uniformly  exponentially  stable. 

Now,  let  us  consider  (ri,r2)  =  (1  -  where  /:  >  0  is  an  odd  integer.  For  k 

large,  this  represents  a  small  perturbation  of  the  delays.  For  this  value  of  the  delays,  it  is 
easy  to  verify  that  the  equation  (1.2)  has  a  solution  Afc  =  i(k  +  ^)Tr.  The  corresponding 
eigenfunction  is  a  solution  of  (1.1)  which  is  periodic  of  period  27r/At.  In  partiuclar,  this 
implies  that  r(cr(5(j__^  2)(1))  ^  1-  There  is  a  general  result  which  we  will  mention  later 
which  shows  that  this  quantity  actually  is  equal  to  1.  For  this  simple  case,  this  fact  also  is 
easy  to  verify  directly.  In  any  case,  we  have  demonstrated  the  following  rather  surprising 
result:  Uniform  exponential  stability  can  be  destroyed  by  a  small  change  in  the  delays. 

It  is  easy  to  modify  the  example  (1.1)  in  sucli  a  way  that,  for  the  delays  (ri ,  r2)  =  (1, 2) 
the  solution  x  =  0  is  uniformly  exponentially  stable  and,  for  the  delays  (ri,r2)  =  (1  — 
2]^, 2),  there  is  a  solution  wliich  becomes  unbounded  at  an  exponential  rate.  It  is  only 
necessary  to  replace  the  coefficients  (A,  i)  by  {\P,  |/?),  where  ^  >  1,  but  sufficiently  close 
to  1. 

We  can  obtain  the  same  type  of  results  for  the  case  when  we  consider  only  small  delays 
since  this  is  a  matter  of  rescaling  time. 

Example  1.2  (Small  delays  in  neutral  equations  can  be  bad).  Consider  the  equation 

(1-3)  -  n)  +  ^a:(<  -  r2)]  =  -yx(t). 

where  7  >  0  is  a  constant.  If  we  choose  initial  data  (p  £  C,  there  is  a  unique  solution 
x(t,<p)  of  (1.3)  which  is  defined  for  all  t  G  IR.  If  we  define  [T(r,,rj)(0v’](^)  = 

Q  €  (— co,0],  then  T^ri,r2)(l'),  t  >  0,  is  a  C°-seinigroup  on  C. 

We  make  use  of  the  following  result  (see  Hale  and  Verduyn-Lunel  (1993)),  which  can 
be  proved  using  the  Laplace  Transform, 

f 

Theorem  1.2. 

(i)  »'(a(r(ri,r,)(l))  =  where  is  deSned  by 
(1-4)  =  ^(ri,r2)  =  sup{ReA  :  =  0}, 
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^7i(»‘i,rs)(-^)  —  -^(l  +  2®  *  2^  *1  + 'y- 

(ii)  re(a(T(r.,r,)(l))  =  r(a(5(r,,r,)(l)).  S(r^,r:,)(t)  is  the  semigroup  for  (1.1). 

Let  us  apply  this  result  for  different  values  of  (rj,  rj).  If  we  choose  the  delays  (rj ,  r2)  = 
(1,2),  then,  from  Example  1.1,  we  know  that  re(a(r(rj,rj)(l))  =  1/v^  <  1-  Also,  it  is 
possible  to  choose  7  so  that  there  is  a  solution  Aq  of  A.y^(r,,r2)(A)  =  0  with  ReAo  > 
— (ln2)/2.  From  Theorem  1.2,  this  means  that  the  exponential  decay  rate  of  the  semigroup 
T(i^2)(t)  is  determined  by  an  element  of  the  point  spectrum  and  not  from  the  essential 
spectrum. 

Now,  let  us  choose  (ri,r2)  =  (1  —  2Fh>^)'  Example  1.1,  we  know  that 

r(<7(Sfi _ s— ,2)(i^))  ~  From  Theorem  1.2,  this  implies  that  re(<7(r(r,,rj)(l))  =  1-  The 

uniform  exponential  stability  of  the  origin  is  destroyed  by  a  small  change  in  (ri,r2)  and  it 
has  nothing  to  do  with  the  constant  7. 

As  in  Example  1.1,  we  can  obtain  exponential  growth  by  changing  slightly  the  coeffi¬ 
cients  in  the  difference  equation. 

We  can  also  achieve  the  same  effect  with  the  delays  being  small  if  we  simply  rescale 
time  t  t/e  and  changing  the  constant  7  i-»  €7,  where  c  is  a  small  parameter.  The 
resulting  equation  for  e  =  0  is  the  ordinary  differential  equation  i  =  —(7/2)1,  which 
is  imiformly  exponentially  stable.  Changes  in  the  delays  as  above  destroy  this  type  of 
stability. 

Why  is  it  to  be  expected  that  such  drastic  changes  in  the  dynamics  as  exhibited  in 
Examples  1.1  and  1.2  can  result  from  small  changes  in  the  delay?  General  resiJts  on 
continuous  dependence  of  the  spectrum  of  boimded  operators  on  a  Banach  space  assume 
that  the  perturbations  are  bounded  .  Even  the  radius  of  the  essential  spectrum  cannot 
increase  drastically  under  such  perturbations.  H  we  consider  the  operator  D(^rx,ri)  '•  C  C 
defined  by 

D^p  =  V5(0)  +  ^<p(-ri )  -1-  ^<,3(-r2), 

then  it  is  not  bounded  as  a  function  of  (rj,r2)- 

We  now  put  the  conclusions  from  these  examples  in  a  more  general  context.  First,  we 
consider  the  system  of  neutral  differential  difference  equations 

(1.5)  ^[x{t)-'^jLiaMt-rj)]  =  i:^=,Qbjx(t-rj),  • 

where  ro  =  1,  0  <  rj  <  co,aj,bj,j  =  1,...,M  are  constants.  Without  loss  of  generality, 
we  can  assume  that  n  <  r2  <  ...  <  rjjf.  Let  r  =  (ri, . . .  ,rAf),  a  =  (ai,...,aM)»  — 

(bo,  61, ... ,  5m),  and  let  Tr(<)  =  T^r,a,b){t)  be  the  C°-semigroup  on  C  defined  by  (1.5).  Also, 
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let  5(r,a)(0  C°-semigroup  on  Co  =  C'o,(r,a)  =  W  ^  C  :  (^(0)  -  E^jaj(^(-rj)  =  0} 

defined  by  the  difference  equation 

(1.6)  x(<)-E;i,G,x(<-r,)  =  0. 

Theorem  1.2  holds  in  this  more  general  situation  (see  Hale  and  Verduyn-Lunel  (1993)). 
Theorem  1.3. 

(i)  r(a(r(r,a,6)(l)))  =  where  ^  is  deGned  by 


(1.7) 


^  =  Pir,a,b)  =  sup{ReA  :  A(^r,c,b)i^)  - 
A(^r,a,b)W  =  <^[1  -  • 

(ii)  re(cr(5(r,a)(l)))  =  e“,  where  ot  is  deGned  by 


(1.8)  Q  =  Q(r,a)  =  sup{  Re  A  :  A(r,a)(A)  =  0}, 

(Hi)  r,(<7(T(.,..„(l)))  =  .■(<-(S(r,.)(l)))- 

From  this  result,  it  is  clear  that  the  eigenvalues  of  the  difference  equation  (1.6);  that 
is,  the  zeros  of  A(r,a)(A)  =  0,  play  a  fimdamental  role  in  the  discussion  of  the  asymp¬ 
totic  behavior  of  the  solutions  of  (1.5)  when  we  subject  the  equation  to  perturbations  in 
the  delays.  We  now  discuss  some  properties  of  these  eigenvalues  and  relate  them  to  the 
semigroup  5(r,a)(l)-  The  presentation  follows  Avellar  and  Hale  (19S0)  and  includes  some 
results  obtained  previously  by  Moreno  (1973),  Henry  (1974),  Melvin  (1974),  Hale  (1975) 
and  Silkowski  (1976). 

Theorem  1.4. 

(i)  r(<7(5(r,a)(l)))  is  continuous  in  a. 

(ii)  r(cr(5(r,a)(l)))  is  lower  semicontinuous  in  r. 

(Hi)  r(cr(5(r,a)(l)))  is  continuous  in  r  if  the  components  of  r  are  rationally  independent. 

(iv)  If  the  components  of  r  are  rationally  independent,  then  r(<7(5(r, «)(!)))  =  c'*®, 
where 

1  = 

(v)  Re<7(5(r,a)(l))  c  [pM,o],  where  a  is  deGned  in  (1.8)  and  pM  satisfies 
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Furthermore,  if  the  components  of  r  are  rationally  independent,  then  \pM,Po]  is  the 
smallest  closed  interval  containing  Recr(S(r,a)(l))- 

Corollary  1,1.  Exponential  stability  of  the  difference  equation  (1.6)  is  preserved  under 
small  pertubations  in  the  delays  if  and  only  if 

(1-9)  Sjlja,l<l. 

The  proof  is  an  easy  consequence  of  Theorem  1.4.  In  fact,  from  Theorem  1.4  (iv), 
if  the  components  of  r  are  rationally  independent,  then  exponential  stability  of  (1.6)  is 
equivalent  to  having  po  <  0.  Also,  from  the  definition  of  po,  we  see  that  po  <  0  if  and  only 
if  (1.9)  is  satisfied.  Since  the  rationals  are  dense  in  R,  we  obtain  the  result. 

Corollary  1.2.  If  exponentially  stability  of  (1.5)  is  to  be  preserved  under  small  pertur¬ 
bations  in  the  delays,  then  it  is  necessary  and  sufficient  that  (1.9)  is  satisfied. 

The  proof  is  a  direct  consequence  of  Theorem  1.3  (i),  (iii),  Corollary  1.1  and  Rouche’s 
Theorem. 

Example  1.1  (Revisited)  In  Example  1.1,  we  considered  in  some  detail  the  behavior  of 
the  zeros  of  the  function 

In  this  case,  it  is  easy  to  see  that  po(r)  =  0  for  all  r.  Also,  if  the  components  of  r  axe 
rationally  independent,  Theorem  1.4  implies  that  [p2(r),0]  is  the  smallest  closed  interval 
containing  Rea(5(r,a)(l)).  A  simple  computation  shows  that  P2(l,2)  =  -ln2.  Since  p2{r) 
is  continuous  in  r,  it  follows  that  the  smallest  closed  interval  containing  Re  (T(S(r  «)(!))  is 
[  In  2,0]  for  any  r  close  to  (1,2)  which  has  rationally  independent  coefficients. 

Example  1.3  (Delays  approaching  zero).  Let  us  consider  the  equation 

h(A,  c,  e)  =  1  -  2ce“'^“  +  =  0, 

where  0  <  Cj  <  cj,  and  investigate  the  behavior  of  the  solutions  as  e  '=  (cj,  €2)  0.  Avellar 

and  Hale  (1980)  prove  the  following  result:  If  -2\c\  <!-<?<  2|c|,  then  p2(e)  -4  -00, 
Po(e)  +00  aj  e  -4  0.  From  Theorem  1.4,  if  the  components  of  c  are  rationally  indepen¬ 
dent,  then  [p2(c),  Po(c)]  is  the  smallest  closed  interval  containing  Re  <r(5(r,o)(l)).  Thus,  if  c 
satisfies  the  above  restrictions  then  the  smallest  closed  interval  containing  Re <t(5(,.  „)(!)) 

approaches  the  interval  (—00,00)  as  € —4  0. 

It  also  is  possible  to  show  that,  if  we  choose  |cl  <  1  and  |c|£2/2ei  =  1,  then  the  closure 
of  the  set  Re<7(5(r,a)(l))  is  equal  to  the  interval  [p2(£),Po(e)]. 

On  the  other  hand,  if  we  choose  £3  =  2€i  and  \c\  <  1,  we  see  that  h(A,c,  c)  =  0  if  and 
only  if  Re  A  =  (l/£i)ln|cl,  which  approaches  -00  as  ej  -4  0.  These  remarks  show  that. 
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if  we  change  the  slope  of  the  hne  along  which  c  — »  0,  then  the  structure  of  Re r,o)(l)) 
changes  drastically. 


In  the  applications,  it  is  not  ailways  the  case  that  the  delays  vaxy  in  an  independent 
manner.  To  be  precise,  we  need  some  notation.  Let  7^-  =  (7ji, . •  •  ,7jm)7  'tjk  nonnegative 
integers  for  all  ;  =  1, . . . , N,  =  1, . . . , M,  70  =  (0, . . .  ,0)  €  (E+)^,  7^  •  r  =  2:^i7jtri, 
a  =  (c], . . .  ,a/y/)  €  6  =  (61,. . .  ,6yv)  €  E^,  and  assume  that 


(1.10)  0  <  7i  •  r  <  72  -  r  <  , . .  <  7W  •  r. 

We  consider  the  NDDE 

(l.U)  ily(<)  +  -  7y  ■  r))  =  -  7,  ■  r). 

The  characteristic  equation  for  (1-11)  is 


(1.12) 

A(r,a,i)('^)  =  AA(r,a)(A)  +  g{X,  6,  r)  =  0, 

where 

(1.13) 

A(r.a)(A)  =  l  +  Ejtjaye-^^^-, 

(1.14) 

g{X,b,r)  =  ^^=,bje-^-^rr_ 

Along  with  (1.12), 

we  consider  also  the  difference  equation 

(1.15) 

y(t)  +  Ejliajy(t  -  7j  •  r)  =  0. 

As  before,  let  T(r,a,6)(<)  be  the  C°-semigToup  on  C  defined  by  (1.11)  and  let  5(r,a)(<) 
be  the  C°-semigroup  on  Co  =  Co,(r.a)  =  {(p  e  C  :  v?(0)  +  •  r^)  =  0}  defined 

by  (1.15). 

For  (1.11)  and  (1.15),  the  Theorem  1.3  is  valid  as  well  as  Theorem  1.4,  Parts  (i),  (ii) 
and  (iii).  The  other  statements  in  Theorem  1.4  are  not  true  unless  N  =  M  since  the  delays 
are  not  varying  in  an  independent  way.  The  computation  of  r((7(5(,.,a)(l)))  in  the  general 
case  is  much  more  difficult.  To  xmderstand  how  to  find  this,  we  introduce  some  additional 
notation. 

Let 

(1.16) 

Z{a,  r)  =  {  Re  A  :  A^r.a){X)  =  0}, 

and  Z(a,r)  =  Cl Z{a,r),  the  closure  of  Z(a,r). 
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We  first  remark  that,  if  the  equation  A(r,a)(M  +  =  0  some  real  numbers  /x,  i/, 

then  we  must  have,  for  0  <  j  <  N, 

(1.17) 

We  define  the  numbers  pj  =  pj(a,r),  0  <j  <  AT,  if  they  exist,  by  the  relations 

(1.18)  = 

where  co  =  1.  It  is  easy  to  verify  that  and  pq  are  uniquely  defined  and  have  the 
property  that  pyv  =  po  if  N  =  1  and  ppj  <  po  if  N  >  2.  Furthermore,  both  /JA((a,r)  and 
Po{o;^)  are  continuous  in  a,r.  Also,  from  (1.17),  it  is  clear  that  we  have  the  following 
relation: 

(1-19)  Z(a,r)  C  [/9N(a,r),po(a,r)]. 

Relation  (1.19)  gives  an  estimate  on  r(<7(S(r,a)(l))),  but  it  may  be  very  inacciarate. 

The  following  result  gives  a  a  more  explicit  characterization  of  the  set  Z(a,r).  The 
proof  is  an  easy  consequence  of  Kronecker’s  Theorem. 

Theorem  1.5.  If  6  =  {6i,...,6m),  0<6j  <  2-k,  ;  =  1,2,...,M,  and 
(1-20)  H{p,e,a,r)  =  1  +  ^ 

and  the  components  of  r  are  rationally  independent,  then  p  €  Z(a,r)  if  and  only  if  there 
is  a  6  such  that  H{p,  9,  a,  r)  =  0. 

Corollary  2.1.  For  any  r  €  (IR^)^,  Z{a,r)  is  the  union  of  a  £nite  number  of  intervals. 

Proof.  K  the  components  of  r  are  rationally  independent,  then  Z(a,  r)  is  characterized  by 
the  solutions  of  H{p,d,a,r)  -  0,  where  H  is  defined  in  (1.20).  Since  these  solutions  are 
einalytic  varieties,  it  is  impossible  to  have  the  following  property:  there  exists  ape  Z{a,r), 
{pj}  e  Z{a,r),  pj  p  as  j  —*  oo,  (pj+i,Pj)  fl  Z(a,r)  =  0.  This  proves  the  result  when 
the  components  of  r  are  rationally  independent. 

For  any  r  €  (IR^)*^,  there  exists  a  ^  e  (IR©  )’  for  some  integer  q  such  that  the 
components  of  ^  are  rationally  independent.  Apply  the  previous  result  to  P  to  complete 
the  proof. 

Example  1.4.  Let  us  consider  the  characteristic  equation 
(1.21)  Ao(p,  h)  =  1  +  fee-''*  - 
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where  a  >  0,  >  0  are  constants.  FVom  Theorem  1.5,  p  €  Z{h)  =  Cl  {Re  /z :  Ao(/z,  h)  =  0} 
if  an  only  if  there  exist  61,62  £  [0, 27r]  such  that 

(1  99)  gZidj  _  2(p-ha)  1  + 

If  we  define 


(1.23) 


h{p)  = 


1  _  e-2(p+a) 

1  +  e-2(p+<i)  ’ 


and,  if  there  is  a  /?  such  that 

=  h{p), 

then  the  right  hand  side  of  (1.22)  has  modulus  1,  which  implies  that  p  £  Z{h).  As  a 
consequence,  for  any  p  >  0,  if  k  >  koip),  then  we  can  find  an  /i  =  h(p,k)  >  0  such  that 
ke  =  kQ(p);  that  is,  we  have  p  £  Z(h)  and  the  radius  of  the  essential  spectrum  of  the 
semigroup  at  <  =  1  is  larger  than  1.  If  /:  >  1,  then  tliis  implies  that  we  can  choose  a  p  as 
large  as  we  want  and  then  choose  an  k  sufficiently  small  so  that  p  £  Z(h).  If  k  >  ko(0), 
then  there  exist  small  p  and  small  k  such  that  p  £  Z{h). 

If  we  take  the  delay  /i  =  0  in  (1.21),  then  the  solutions  of  (1.21)  must  have  negative 
real  parts  for  any  A:  >  0;  that  is,  vve  have  expontential  stabihty  for  the  semigroup  of  the 
corresponding  difference  equation.  On  the  other  hand,  if  k  1,  we  can  choose  a  p  as  large 
as  we  want  and  then  choose  an  h  sufficiently  small  so  that  p  £  Z{h)-,  that  is,  there  is  a 
solution  of  the  difference  equation  becoming  unbounded  at  a  very  large  exponential  rate. 

In  the  next  section,  we  show'  that  this  example  occurs  also  in  the  characteristic  equa¬ 
tion  arising  from  the  control  of  a  wave  equation  with  delayed  boundary  control. 

For  other  examples,  see  Avellar  and  Hale  (19S0),  where  they  also  discuss  similar 
problems  for  matrix  equations.  Further  discussion  is  contained  in  Hale  and  Verduyn-Lunel 
(1993). 


2.  Delayed  Boundary  Control.  In  recent  years,  there  has  been  considerable  effort 
devoted  to  the  problem  of  stabilization  and  control  of  PDE  through  the  application  offerees 
on  the  boundary.  The  mathematical  theory  is  very  complete  in  the  situation  when  the 
boundary  forces  are  applied  with  no  delays  in  time  (see,  for  example.  Lions  (1988)).  It  has 
been  pointed  out  recently  that  a  small  time  delay  in  the  application  of  the  boundary  forces 
can  lead  to  a  destabilization  of  the  system  (see  Datko,  Lagnese  and  Polis  (1986),  Datko 
(1988),  (1991),  (1994),  Desch  and  Wheeler  (1989),  Hannsgen,  Renardy  and  Wheeler  (1988), 
Logemann,  Robarber  and  Weiss  (1993),  and  the  references  therein).  In  this  section,  we  give 
such  an  example  for  the  wave  equation,  showing  that  the  drastic  changes  in  the  stability 
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properties  occur  essenticJly  for  the  same  reason  that  we  have  noticed  in  the  previous  section 
for  difference  equations. 

Consider  the  linear  wave  equation 

(21)  v}it  +  2axvt-Wj:i  +  bw  =  0,  0  <  z  <  1,  <>  0, 

with  the  boundary  conditions 

(2-2)  u)(0,  t)  =  0,  t)  =  -kwt{l, t  -  h), 

where  n  ^  0,  /:  >  0,  /i  >  0, 6  are  constants.  System  (2.1),  (2.2)  corresponds  to  a  boundary 
stabilization  problem  where  the  control  function  is  kwt{l,t-  h).  It  is  not  difficult  to  show 
that  this  system  generates  a  C°-semigroup  Sh,k{i)  on  the  space  H^(0,  1)xL2(0,  1).  We  want 
to  determine  how  rt(Sh,k{l))  depends  upon  k,h.  If  we  replace  w{x,t)  by  c~“‘u;(z,t)  and 
use  the  fact  that,  for  any  constant  c,  a  term  cw  corresponds  to  a  compact  perturbation  of 
the  differential  operator  then  we  see  that  re(SA.fc(l))  =  where  Sk,k{t) 

is  the  semigroup  generated  by  the  equation 

(2-3)  wtt  -  Wxx  =  0,  0  <  z  <  1,  t  >  0, 

with  the  boundary  conditions 

(2.4)  u>(0,0  =  0,  U’x(l,t)  = -A:e“'‘[-aii)(l,t-/i) +  10^(1,  t-h)]. 

These  equations  can  be  written  in  an  equivalent  form 

(2-5)  ut  =  -Ux,  vt  =  -Ux, 

with  the  boundary  conditions 

(2.6)  .  t;(0,t)=0,  «t(l,t)  = -/:e“'*[-au(l,t-h)  +  t;t(l,t-h)]. 

Let  us  show  that  (2.5),  (2.6)  is  equivalent  to  a  NDDE  with  three  delays.  The  general 
solution  of  the  partial  differential  equation  is 

v{x,t)  -  <p(x  -  <)  +  rl>(x  +  t), 

u(z,  t)  =  <p(x  -t)-  rlf(x  +  <). 

This  implies  that 

2(p(x  -t)  =  v(z,t)  +  u(z,  t) 

2ip{x  + 1)  =  v(i,t)  -  u(z,t). 
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From  these  expressions,  we  deduce  that 

2ip{-t)  =  v{l,t  +  \)  +  u{l,t  +  l) 

2T/>(0  =  v(l,t-l)-u(l,t-l)'. 

Using  (2.7)  and  the  first  boundary  condition  at  t  —  1,  we  deduce  that 

u(l,t)  —  u(l,t  —  2)  =  — v(l,t)  —  v(l,<  —  2). 

Differentiating  this  expression  with  respect  to  t,  using  the  second  boundary  condition  and 
letting  y{t)  =  v(l,t),  we  conclude  that 

^[y{t)  +  -  h)  +  y(t  -  2)  -  ke^^y^t  -2-h)]  =  ake'^^yit  -  h)  -  ake^*^y{t  -2-h). 

at 

From  this  equation,  we  see  that  the  essential  spectral  radius  re(5A,it(l))  of  Sh,k{^)  is 
determined  by  the  supremum  of  the  real  parts  of  the  solutions  of  the  equation 

1  +  =  0. 

From  our  derivations  above,  re(5/i,fc(l))  of  S/i,t(l),  is  determined  by  the  supremum  of  the 
real  parts  of  the  solutions  of  the  equation  (let  X  —  a  =  fi) 

(2.8)  1  +  =  0. 

If  this  supremum  is  negative  (positive),  then  solutions  approach  zero  exponentially  (there 
are  solutions  which  are  exponentially  unbounded)  as  t  —*  oo. 

Equation  (2.8)  is  the  same  as  Equation  (1.21)  in  Example  1.4.  As  a  consequence  of 
our  discussion  in  that  example,  we  conclude  that 

(i)  h  =  0,  /:  >  0  implies  exponential  approzw:h  to  zero. 

(ii)  h>0,0<k<{l  —  e~^“)/(l  +  e“^“)  implies  exponential  approach  to  zero. 

(iii)  fe  >  (1  —  e“^“)/(l  +  c~^“)  implies  that  there  is  a  dense  set  of  /i  >  0  such  that 
there  are  solutions  which  are  exponentially  imbounded. 

(iv) H  ifc  >  1,  then  we  can  choose  a  p  as  large  as  we  want  and  then  choose  an  h 
sufficiently  small  so  that  there  is  a  solution  of  (2.1),  (2.2)  which  is  becomes  unbounded  at 
the  rate  e'’*  as  t  — »  oo. 

As  a  consequence  of  this  observation,  it  follows  that  the  system  could  have  been 
stabilized  with  a  control  which  involves  no  delay  and  then  there  are  arbitrarily  small 
delays  in  the  control  which  lead  to  destabilization.  These  remarks  are  contained  in  Datko 
(1991),  (1994),  but  his  proof  is  somewhat  different. 
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#  .  4 


The  fact  that  the  above  linear  wave  equation  is  equivalent  to  a  NDDE  is  a  consequence 

of  a  much  more  general  theory  (see,  for  example,  Hale  and  Verduyn-Lunel  (1993)  and  the 
references  therein). 
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